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Abstract 


Digital  simulations  were  developed  to  implement  a pitch 
rate  control  system  for  the  F-16  aircraft  engaged  in  aerial 
gunnery.  First,  the  EASY  Modelling  and  Analysis  Program 
by  Boeing  Computer  Services  was  adapted  to  implement  a 
longitudinal  axis  F-16  aircraft,  flight  control  system,  and 
pilot  model.  Comparison  of  closed  loop  system  responses 
indicated  a proposed  pitch  rate  flight  control  configuration 
would  improve  target  tracking  performance.  The  Terminal 
Aerial  Weapon  Delivery  Simulation  (TAWDS)  program  by  McDon- 
nell Douglas  Corporation  was  adapted  for  the  F-16  aircraft. 

A non-linear,  six-degree-of- freedom  aircraft  model,  multi- 
axis flight  control  system,  and  multi-axis  pilot  model  were 
developed  to  demonstrate  target  tracking  capabilities. 

Eight  different  air-to-air  scenarios  were  developed  to  simu- 
late evasive  encounters  with  an  F-4  target  aircraft.  Time 
history  target  tracking  errors  indicated  the  improved 
tracking  performance  of  the  proposed  pitch  rate  flight 
control  configuration  over  the  present  normal  acceleration 
configuration  of  the  F-16  aircraft. 


! 


c 


xiv 


PITCH  RATE  FLIGHT  CONTROL 


FOR  THE  F-16  AIRCRAFT  TO 


IMPROVE  AIR-TO-AIR  COMBAT 


I.  Introduction 


One  of  the  most  challenging  tasks  facing  today's  tacti- 
cal fighter  pilot  is  that  of  air-to-air  gunnery.  When 
engaged  in  aerial  combat  with  an  enemy  aircraft,  today's 
fighter  pilot  must  maintain  an  offensive  role  by  tracking 
his  target.  To  be  successful,  he  must  achieve  a target 
tracking  solution  that  allows  him  to  deliver  his  weapons 
quickly  and  accurately.  All  too  often,  this  task  requires 
more  skill  and  control  precision  than  the  pilot  is  able  to 
provide. 

During  recent  years,  considerable  U5AF  and  industrial 
efforts  have  been  directed  to  developing  advanced  tactical 
aircraft  flight  control  systems  to  improve  weapon  delivery 
accuracy.  The  most  promising  of  these  engineering  efforts 
involves  integration  of  aircraft  flight  and  fire  control 
systems.  The  benefits  of  automatic  flight  control,  coupled 
with  automatic  weapon  delivery,  will  allow  a fighter  pilot, 
while  engaged  in  air-to-air  aerial  combat,  to  select  a 
degree  of  automation  to  assist  him  in  both  flying  his 
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aircraft  and  firing  his  weapons  more  effectively.  This 
could  range  from  a manual  system  to  a fully  automatic  mode 
of  operation.  Preliminary  investigations  have  shown  that 
weapon  delivery  accuracy  can  be  improved  by  coupling  air- 
craft flight  control  and  weapon  delivery  fire  control 
systems  to  relieve  the  fighter  pilot  of  his  ever  increas- 
ing workload  (Ref  1).  However,  the  possibility  also  exists 
of  improving  weapon  delivery  effectiveness  by  conditioning 
flight  control  systems  without  removing  the  pilot  from  his 
primary  tasks.  Manual  flight  control  investigations  of  im- 
proving aerial  gunnery  will  be  considered  in  this  study. 

Flight  Control  System  Background 

The  performance  of  air  superiority  aircraft  such  as  the 
McDonnell  Douglas  F-15  or  the  General  Dynamics  F-16  in  air 
combat  maneuvers  places  unusual  and  heavy  demands  on  the 
flight  control  system.  This  is  true  because  today's  high 
performance  aircraft  are  operated  over  extremely  wide  flight 
envelopes.  In  addition  to  using  the  total  altitude  and  Mach 
range,  the  pilot  exercises  the  aircraft  through  its  full 
angle  of  attack  capabilities  during  air  combat  (Ref  2). 

To  aid  the  pilot  in  his  primary  task,  handling  quality 
specifications  have  been  designed  using  a weighted  combina- 
tion of  pitch  rate,  normal  acceleration,  and  pitching 
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acceleration  criteria  (Ref  3).  Aircraft  flight  test  perform- 
ance ratings  have  indicated  a pilot  preference  of  this 
blended  system  for  normal  cruise  maneuvers.  As  a conse- 
quence, several  systems  have  been  built  which  combine 
pitch  rate  and  normal  acceleration  as  the  feedback  variables. 

One  method  for  mechanizing  the  flight  control  system 
feel/response  is  the  C*  approach  that  was  first  proposed 
by  Boeing  aircraft  design  engineers  (Ref  3).  This  approach 
uses  a linear  blend  of  normal  acceleration,  pitch  rate,  and 
pitch  acceleration.  The  weighted  control  combination  can 
be  described  as  follows: 

C*  = k,  A + k09  + k '0  (1) 

in  z o 

where 

A = normal  acceleration  at  c.g. 
n * 

0 = pitch  rate 
8 = pitch  acceleration 

The  C*  equation  can  also  be  defined  in  g's  where  the 
units  of  kg  are  equivalent  to  a velocity  divided  by  gravity 
(g)  and  k^  is  equivalent  to  the  distance  between  the  linear 
accelerometer  and  the  center  of  gravity  of  the  aircraft 
divided  by  g.  Using  k^  = 1,  Equation  (1)  can  be  written  as: 

3 


where 


U 6 L8 

C*  = A + -£2-  + — 
n g g 


U = the  cross-over  velocity  (approximately  400  ft/sec) 
co 

L = distance  between  linear  accelerometer  and  the 
center  of  gravity  of  the  aircraft. 

Selection  of  the  cross-over  velocity  specifies  the 
operating  point  at  which  the  control  contributions  of  pitch 
rate  and  normal  acceleration  feedback  are  equal.  At  lower 
airspeeds,  such  as  in  landing  approaches  where  the  control 
surfaces  are  relatively  less  responsive  than  at  cruise 
flight,  the  pitch  rate  feedback  is  predominant.  At  air- 
speeds above  the  cross-over  velocity,  in  flight  regions 
where  the  aircraft  control  surfaces  are  relatively  more 
responsive  than  at  lower  airspeeds,  the  normal  acceleration 
feedback  is  predominant.  The  C*  approach  is  convenient  for 
the  mechanization  of  a feel  system  because  of  the  ease  by 
which  the  pitch  rate,  pitch  acceleration,  and  normal  accel- 
erations can  be  measured  (Ref  3). 


Gunsiqht  Technology 

Although  present  technology  has  developed  very  advanced 
flight  control  systems  that  could  be  used  in  the  integration 
of  flight  and  fire  control  systems,  a key  limiting  factor  in 
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advancing  air-to-air  combat  is  target  tracking  systems.  The 
aircraft  gunsight  is  mechanized  to  compute  and  display  to 
the  pilot  the  lead  angle  required  to  hit  the  target.  This 
is  generally  done  by  displacing  an  aiming  reticle  for  the 
required  lead  angle  from  some  gunsight  reference  line  to  the 
target.  Essentially,  if  the  pilot  flies  his  aircraft  so  as 
to  keep  the  reticle  on  the  target,  he  then  is  maintaining 
the  proper  lead  angle  to  achieve  a target  hit  (Ref  2). 

For  many  years,  target  tracking  has  been  accomplished 
using  a disturbed  reticle  sight.  The  most  popular  of  these 
sights  and  the  one  that  is  used  on  most  present  day  fighter 
aircraft  is  the  Lead  Computing  Optical  Sight  System  (LCOSS). 
The  LCOSS  generates  a gunsight  lead  angle  based  on  the 
attacker  aircraft's  own  dynamics.  The  attacker  aircraft's 
own  body  rates,  load  factor,  angle  of  attack,  and  airspeed 
are  used  to  determine  a gunsight  lead  angle  solution.  To 
achieve  a continuous  target  tracking  solution  with  the  LCOSS 
sight,  the  pilot  must  fly  his  aircraft  to  remain  in  the  same 
plane  of  motion  as  the  target  aircraft. 

In  contrast  to  the  LCOSS,  a new  director  gunsight  is 
presently  being  developed  and  demonstration  flight  test 
programs  are  scheduled  to  begin  in  the  near  future.  Gen- 
erally, the  director  sight  incorporates  actual  measurements 


of  target  aircraft  motion  to  determine  a gunsight  solution. 
Line  of  sight  angle  rate  and  position  measurements  are  used 
in  the  director  sight  instead  of  own-ship  body  rates  as  in 
the  LCOSS.  The  director  system  employs  a tracking  device 
such  as  an  angle  tracking  radar  or  an  electro-optical  tracker 
to  measure  target  motion.  With  actual  target  measurements, 
the  director  system  incorporates  a Kalman  filter  to  determine 
the  expected  future  target  position  (Ref  4). 


Objective 

The  primary  objective  of  this  study  will  be  to  compare 
the  target  tracking  performance  of  a director  gunsight  im- 
plementation for  manual  flight  control  involving  two  air- 
craft flight  control  configurations.  Present  handling  qual- 
ities specifications,  supported  by  Cooper-Harper  pilot  rat- 
ings have  indicated  that  normal  acceleration  feedback  may 
be  effective  for  most  flight  conditions.  However,  a pitch 
rate  flight  control  scheme  will  be  investigated  to  improve 
aerial  gunnery  during  air-to-air  combat.  Since  the  director 
gunsight  incorporates  actual  target  measurements  of  angular 
error  and  angular  error  rate,  it  seems  appropriate  that  a 
pitch  rate  control  feedback  scheme  could  be  employed  to 
provide  improvements  in  manual  target  tracking  systems. 
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Plan  of  Attack 

4 * 

To  begin  the  investigation  of  pitch  rate  control,  it 
was  necessary  to  select  an  aircraft  whose  flight  control 
system  incorporates  the  C*  concept.  Selection  of  the  F-16 
flight  control  system  as  a candidate  system  is  discussed  in 
Chapter  II.  Considerations  of  a flight  condition  and  air- 
craft  characteristics  are  also  included  in  Chapter  II. 

Implementation  of  a representative  F-16  pilot  model  is 
discussed  in  Chapter  III.  This  analytical  representation 
allowed  a closed  loop  system  to  be  established  for  man-in- 
the-loop  simulations  necessary  for  manual  flight  control 
evaluation. 

A digital  simulation  of  the  F-16  aircraft  dynamics, 
flight  control  system,  and  pilot  model  is  developed  in 
Chapter  IV.  Analysis  of  the  present  F-16  flight  control 
system  is  described  along  with  the  investigation  of  a 
predominately  pitch  rate  control  configuration. 

A second  digital  simulation  program  is  discussed  in 
Chapter  V.  A non-linear  six-degree-of-f reedom  aircraft, 
flight  control  model,  and  multi-axis  pilot  model  is  de- 
veloped to  provide  a closed  loop  simulation.  To  provide 
realistic  target  tracking  encounters,  a series  of  eight 


different  air-to-air  combat  scenarios  is  developed. 


The  results  of  the  air-to-air  encounters  using  a direc- 
tor gunsight  implementation  is  included  in  Chapter  V.  Time 


history  simulation  data  is  generated  to  measure  the  target 
tracking  performance  of  both  the  present  normal  acceleration 
flight  control  configuration  of  the  F-16  aircraft  and  the 
new  proposed  pitch  rate  flight  control  system. 
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I I . Aircraf t Selection 


Select!  in  Criter ia 

The  F-lo  fighter  aircraft  built  by  General  Dynamics 
Corporation  was  selected  as  the  baseline  aircraft  simula- 
tion model.  The  F-16  was  chosen  because  it  represents  the 
present  state  of  the  art  in  fighter  aircraft  design.  Its 
fly-by-wire  flight  control  system  enabled  design  engineers 
to  harness  a basically  unstable  aircraft  and  obtain  unpre- 
cedented flight  performance.  The  design  of  this  flight 
control  system  incorporates  the  C*  concept  discussed  in 
Chapter  I (Ref  3).  The  F-16  flight  control  system  incorpo- 
rates angle  of  attack,  pitch  rate,  and  normal  acceleration 
feedback.  As  the  C*  concept  implies,  normal  acceleration 
feedback  is  predominant  at  cruise  airspeeds.  A blending  of 
normal  acceleration  and  pitch  rate  feedback  is  employed  in 
the  longitudinal  axis  control  system.  In  addition,  angle  of 
attack  (i.e,  alpha) is  fed  back  to  the  flight  control  system 
to  aid  stability  and  achieve  alpha  limiting  at  high  angles 
of  attack.  This  unique  configuration  makes  the  F-16  aircraft 
a very  likely  candidate  to  examine  various  flight  control 
configurations  and  incorporate  these  into  a manual  flight /fire 
control  system  evaluation  (Ref  5). 
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In  addition  to  this  interesting  flight  control  con- 
figuration, the  F-16  was  selected  because  of  the  ongoing 
joint  programs  between  the  Air  Force  Flight  Dynamics  Labora- 
tory and  General  Electric  Company.  Their  continuing  investi- 
gations of  integrating  flight  and  fire  control  systems  has 
recently  included  a manned  simulation  program  using  the  F-16 
as  a baseline  aircraft  model  (Ref  1). 

Aircraft  Model  Description 

To  complement  the  efforts  of  the  Air  Force  Flight 
Dynamics  Laboratory  and  General  Electric  Company,  realistic 
scenarios  were  developed  for  the  simulation.  The  air-to- 
air  encounters  were  set  with  the  F-16  aircraft  in  its  clean 
configuration  at  a cruise  airspeed  of  .8  Mach  and  altitude 
of  20,000  feet  M3L.  It  was  from  this  cruising  flight  con- 
dition that  the  attacker  aircraft  would  engage  the  enemy 
target.  Assuming  that  the  pilot  had  fired  his  two  available 
Sidewinder  missiles,  he  was  equipped  with  only  his  20-milli- 
meter M-61  conventional  cannon  with  which  to  continue  the 
engagement . 

To  validate  the  aircraft  dynamics  of  the  modelling 
programs,  F-16  aircraft  dynamic  simulation  data  obtained 
from  the  Air  Force  Flight  Dynamics  Laboratory  LAMARS  faci- 
lity was  selected  as  the  desired  test  case.  A digital 


program  by  John  Griffin  (Ref  6)  provided  aerodynamic  data 
for  computer  validation  of  selected  aircraft  configurations 
and  flight  conditions  for  the  F-16  manned  simulations  in 
the  LAMARS  facility.  This  data  serves  as  a basis  for  the 
F-16  digital  simulations  to  be  developed.  Table  I des- 
cribes the  F-16  aircraft  model  characteristics  of  the  test 
case  selected.  A more  complete  aircraft  model  description 
and  detailed  listing  of  the  flight  condition  stability 
derivatives  are  shown  in  Appendix  B. 
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Table  I 

F-16  Aircraft  Model  Characteristics 
(Ref  G) 


Fliqht  Condition 

Altitude 

Airspeed 

Dynamic  pressure 

Air  Density 

20, 000  ft 
.8  Mach 

(829.5  ft/sec  0 20,000  ft) 
436.06  lbs/ ft  „ 

.001267  slug- ft 

Aircraft 

Clean  configuration 

Gross  weight 

19,000  lbs 

Mass 

590.5  slugs 

c.g.  location 

33.92%  moan  aerodynamic 

chord  (MAC) 

Chordlength 

11.32  ft 

Wing  span 

30  ft  9 

Wing  area 

300  ft 

Moments  of  Inertia 

XX- ax is 

9007.5  slugs- ft0' 

YY-axis 

49956.  slugs- ft 

ZZ-axis 

56770.  slugs-ft“ 

XZ-axis 

198.0  slugs-ft" 

Trim  Fliqht  Condition  Parameters 

Load  factor 

1 g (32.174  ft/sec2) 

Flight  path  angle 

0 degrees 

Anglo  of  attack 

2.1039  degrees 

Stabilator 

-1.9123  degrees 

( elevator ) 

(up  deflection) 

Armament 

M-61  gun 

20  mm  ammunition 

- 

3400  ft /sec  muzzle 

velocity 

* 

6000  rounds /minute 
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Basic  F-4  Aircraft  Pilot  Model  Description 

In  order  to  implement  a closed  loop  system  performance 
analysis,  an  F-lo  multi-axis  pilot  model  was  required. 

Manned  simulation  efforts  by  McDonnell  Douglas  Corporation 
have  been  instrumental  in  the  development  of  an  analytical 
pilot  model  for  the  F-4E  aircraft.  Their  mathematical  model 
was  developed  from  target  tracking  data  produced  by  two 
DSAF  pilots  flying  in  the  MCAIR  flight  simulator  (Ref  7).  In 
simulating  air-to-air  weapon  delivery,  a data  base  was  provi- 
ded by  measuring  aircraft  tracking  time  histories,  pilot 
tracking  task  responses,  statistical  tracking  performance 
and  weapon  delivery  performance.  The  pilots  were  required 
to  track  a target  aircraft  through  programmed  maneuvers 
while  their  tracking  performance  responses  were  documented. 

Although  efforts  of  McDonnell  Douglas  were  directed 
to  developing  an  F-4E  aircraft  pilot  model,  they  determined 
from  time  history  data  that  pilot  elevation  and  traverse 
tracking  error  characteristics  were  similar  regardless  of 
the  pilot,  the  weapon  delivery  task,  aircraft  flying  quali- 
ties, or  sight  system  characteristics.  The  results  of  the 
McDonnell  Douglas  study  indicated  that  elevation  tracking 
error  contained  two  predominant  modal  components.  This  was 
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due  to  the  pilot's  interaction  with  the  aircraft's  short 
period  dynamics  and  the  sight  dynamics.  Both  frequency 
components  were  observed  to  exhibit  a limited  or  lightly 
damped  response.  Their  study  of  pilot  responses  indicates 
that  the  longitudinal  pilot  model  treats  the  pilot  as  a 
proportional-plus-derivative  observer  of  the  tracking  error. 
The  pilot  threshold  limit  of  error  rate  is  indicated  by  use 
of  a dead  zone  in  the  error  rate  channel.  This  proportional- 
plus-derivative  observer  of  elevation  angular  error  results 
in  a tracking  error  projected  a time  interval  into  the 
future.  This  projected  error  is  then  coupled  with  the  output 
of  a low  pass  filter  to  determine  the  pilot's  rate  input 
to  the  control  stick.  Integrating  this  control  stick  rate 
provides  the  control  stick  position  or  stick  force  that 
determines  the  pilot  model  response  to  the  flight  control 
system.  The  pilot  model  schematic  is  shown  in  Figure  1. 

Just  as  in  the  longitudinal  pilot  model,  the  lateral- 
directional  model  is  based  on  the  assumption  that  the  pilot 
acts  as  a proportional-plus-derivative  observer  of  the  tra- 
verse tracking  error  with  the  dead  zone  on  the  error  rate. 
However,  additional  visual  ques  that  the  pilot  may  perceive 
in  traverse  tracking  are  incorporated  in  the  lateral-direc- 
tional pilot  model.  This  includes  feedback  of  attacker 


Figure  1.  Multi-Axis  Pilot  Model  for  Air-to-Air  Aerial  Gunnery 
(Longitudinal  and  Lateral-Directional  Axes)  (Ref  7) 


aircraft  bank  angle  relative  to  the  target  aircraft.  This 
additional  input,  which  incorporates  a dead  band  for  thres- 
hold limits,  is  used  to  modify  the  lateral  stick  commands 
based  on  traverse  tracking  error. 

The  schematic  diagram  of  the  longitudinal  and  lateral- 
directional  multi-axis  pilot  model  developed  by  McDonnell 
Douglas  is  shown  in  Figure  1.  A root  locus  stability 
analysis  of  the  pilot  model  transfer  function  was  used  to 
determine  how  the  parameters  of  the  pilot  model  affect 
closed  loop  system  stability  and  control  during  weapon 
delivery.  Both  the  director  sight  and  the  lead  computing 
optical  sight  system  were  used  in  the  air-to-air  gunnery 
investigation  and  pilot  model  validation.  Tracking  perform- 
ance time  histories  and  frequency  responses  of  the  pilot 
model  performing  weapon  delivery  tasks  closely  represented 
those  of  the  human  pilot  in  the  manned  simulations  (Ref  7). 

F-16  Aircraft  Pilot  Model  Description 

Unfortunately,  the  gain  parameters  of  the  pilot  model 
developed  by  McDonnell  Douglas  for  the  F-4  aircraft  configu- 
ration was  not  sufficient  for  an  F-16  simulation.  Since 
this  basic  model  was  not  compatible  with  the  F-16  aerodynamic 
and  flight  control  characteristics,  an  effort  was  made  by 
General  Dynamics  to  adapt  this  basic  F-4  pilot  model  to  the 
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aircraft  bank  angle  relative  to  the  target  aircraft.  This 
additional  input,  which  incorporates  a dead  band  for  thres- 
hold limits,  is  used  to  modify  the  lateral  stick  commands 
based  on  traverse  tracking  error. 

The  schematic  diagram  of  the  longitudinal  and  lateral- 
directional  multi-axis  pilot  model  developed  by  McDonnell 
Douglas  is  shown  in  Figure  1.  A root  locus  stability 
analysis  of  the  pilot  model  transfer  function  was  used  to 
determine  how  the  parameters  of  the  pilot  model  affect 
closed  loop  system  stability  and  control  during  weapon 
delivery.  Both  the  director  sight  and  the  lead  computing 
optical  sight  system  were  used  in  the  air-to-air  gunnery 
investigation  and  pilot  model  validation.  Tracking  perform- 
ance time  histories  and  frequency  responses  of  the  pilot 
model  performing  weapon  delivery  tasks  closely  represented 
those  of  the  human  pilot  in  the  manned  simulations  (Ref  7). 

F-16  Aircraf t Pilot  Model  Descr iption 

Unfortunately,  the  gain  parameters  of  the  pilot  model 
developed  by  McDonnell  Douglas  for  the  F-4  aircraft  configu- 
ration was  not  sufficient  for  an  F-16  simulation.  Since 
this  basic  model  was  not  compatible  with  the  F-16  aerodynamic 
and  flight  control  characteristics,  an  effort  was  made  by 
General  Dynamics  to  adapt  this  basic  F-4  pilot  model  to  the 
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F-16  aerodynamic  and  flight  control  characteristics.  Their 
approach  was  first  a pilot  model  gain  study  followed  by  a 
comparison  of  flight  test  data  with  manned  simulation  re- 
sults. Their  study  provided  a multi-axis  pilot  model  that 


produced  fairly  adequate  tracking  performance  for  a stabi- 
lized 5 g target  aircraft  model.  Table  II  indicates  the 
parameter  gain  values  of  the  multi-axis  pilot  model  adapted 
for  this  simulation.  The  parameters  of  the  pilot  model 
shown  are  for  employing  the  director  sight  system.  Pilot 
performance  variations  while  employing  the  LCOS  system  are 
included  by  the  gain  changes  as  indicated  by  the  asterisk 
(Ref  8). 


Table  II 

Parameter  Values  of  the  F-15  Pilot  Model 
Implementing  Director  Sight  (Ref  8) 


r — — — — — i 

T to) 

.05  1.0 

C 

0.6 

*S  *L  kp* 

.125  .125  .25 

K 

♦ 

.0573 

« 

kail  krtr 

.10  1.5 

*S?UD 

0.0 

DB-1  DB-2 

lmr/sec  5 deg 

DB-3 

2.5mr /sec 

#To  implement  the 

LCOS, 

the  above  parameters 

remain  the  same 

except : 

KAIL  ■ -17 
% * -20 

• 
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IV.  Development  of  the  Analysis  Model 


EASY  Program  Discuss  ion 

To  aid  in  the  analysis  of  the  present  flight  control 
configuration  and  implementation  of  a pitch  rate  controller, 
the  EASY  Modelling  and  Analysis  Program  by  Boeing  Computer 
Service  was  adapted  (Ref  9).  The  EASY  program  package  con- 
sists of  two  programs  which  allow  the  modelling  and  analy- 
sis of  dynamic  aircraft  systems  in  both  steady  state  and 
dynamic  behavior.  The  first  of  these  is  the  EASY  Model 
Generation  program.  This  pre-compiler  program  accepts  model 
description  instructions  from  the  programmer  and  from  these 
instructions  generates  a FORTRAN  model  of  the  aircraft  system. 
The  output  of  the  EASY  model  program  is  a complete  system 
model  description  including  a computer  generated  schematic 
diagram  showing  inter-connections  between  the  components  of 
the  constructed  model.  Standard  EASY  components  used  in- 
clude aircraft  modelling  components  and  control  system 
components.  The  computerized  model  was  analyzed  using  the 
linear,  non-linear,  steady  state,  and  dynamic  techniques 
available  in  the  EASY  Analysis  program.  FORTRAN  statements, 
referred  to  as  "program  commands",  allowed  the  analysis  to 
include  non-linear  simulation,  steady  state  calculations, 
linear  model  generation  from  the  original  non-linear  model, 
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eigenvalue  calculation,  root  locus  analysis,  transfer 
function  calculation  and  several  other  dynamic  analysis 
techniques . 

F-16  EASY  Model 

The  EASY  Model  Generation  program  was  begun  by  develop- 
ing a schematic  diagram  of  the  longitudinal  axis  of  the  F-16 
flight  control  system.  The  fold-out  of  Figure  2 indicates 
the  present  F-16  aircraft  control  configuration.  Simplifi- 
cations of  the  pitch  axis  system  (upper  left  of  Figure  2) 
were  made  for  the  EASY  model  program.  The  aircraft  flight 
condition  chosen  was  .8  Mach  at  20,000  feet  and  the  aircraft 
model  configuration  was  that  for  cruise  flight  with  no 
pitch  trim  nor  pitch  autopilot  included.  Trailing  edge 
flaps  were  not  implemented  and,  because  very  little  control 
blending  occurs  at  cruise  airspeeds,  the  differential  tail 
deflection  of  the  F-16  aircraft  was  also  not  modelled. 
Instead,  it  was  assumed  that  the  aircraft  exhibits  conven- 
tional elevator  deflections.  Assuming  rigid  body  dynamics, 
the  high  frequency  structural  filters  were  also  omitted. 

The  resulting  F-16  longitudinal  flight  control  system 
modelled  for  the  EASY  analysis  program  is  shown  in  the 
schematic  diagram  of  Figure  3 (Ref  10). 
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LIGHT  CONTROL  SYSTEM  FUNCTIONAL  BLOCK  DIAGRAM 


Figure  2.  Present  F-16  Flight  Control  System  Diagram  (Ref  11) 
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Only  standard  components  of  the  EASY  program  were 
necessary  to  complete  the  flight  control  system  modelling. 

The  dash  boxes  around  the  elements  in  the  control  system 
diagram  of  Figure  3 indicate  each  standard  component  used. 

For  example,  the  FU  components  indicate  table  look-up 
functions.  These  include  both  the  pitch  command  gradient 
and  gain  scheduled  variables  that  are  functions  of  dynamic 
and  static  pressure.  The  LA  and  LG  components  were  used  to 
model  first  order  lag  transfer  functions.  First  order  lead- 
lag  transfer  functions  wore  modelled  with  the  LE  component. 
The  multiply  and  add  components,  MA  and  MC,  were  used  to 
model  the  summing  junctions.  Second  order  transfer  functions 
such  as  in  the  elevator  actuator  were  modelled  using  the 
TF  component.  Limiting  functions  or  saturation  function 
components,  SA,  were  used  to  regulate  the  pilot  commanded 
maximum  g forces,  angle  of  attack,  actuator  rates,  and 
elevator  deflection. 

In  addition  to  specifying  the  F-1G  flight  control 
system  components,  the  aircraft  dynamic  modelling  was 
necessary.  Description  of  the  aircraft  motion  centered 
around  the  standard  components  AV,  LO,  and  SD.  The  AV 
component  uses  the  aircraft  states  to  compute  aerodynamic 
variables  such  as  angle  of  attack,  airspeed,  body  rates,  etc. 
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The  longitudinal  aerodynamic  component,  LO,  computes  the 
longitudinal  aerodynamic  forces  and  moments.  The  six-degree- 
of-freedom  equations  of  motion  component,  SD,  contains  the 
rigid  body  dynamic  equations  for  integrating  the  aircraft 
states  and  is  driven  by  the  aerodynamic  variables  generated 
in  LO. 

In  addition  to  the  aircraft  and  flight  control  system 
model,  the  longitudinal  pilot  model  described  in  Chapter  III 
was  also  implemented  into  the  EASY  program.  Again,  standard 
components  were  used  to  complete  the  pilot  model  description 
as  shown  in  the  schematic  diagram  of  Figure  4.  This  pilot 
model  implementation  incorporates  parameter  requirements 
for  the  director  sight.  Although  the  description  appears 
different  from  that  of  Figure  1,  p.  16,  a mathematically 
equivalent  model  is  shown.  The  time  delay  component  of 
Figure  1 has  been  incorporated  in  both  the  proportional  and 
derivative  channels  of  the  longitudinal  pilot  model.  In 
addition.  Figure  4 indicates  the  pseudo  target  tracking  task 
of  the  pilot  model  for  the  EASY  analysis.  A closed  loop 
system  was  achieved  by  feeding  back  the  attacker  aircraft 
pitch  angle.  This  pitch  angle  was  compared  to  a reference 
angle  to  allow  performance  evaluations  of  system  response 
to  step  inputs.  The  aircraft  pitch  angle  was  chosen  since 
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tudinal  Pilot  Model  Schematic 


it  is  angular  target  measure  that  is  observed  by  the  pilot 
in  his  target  tracking  task.  Unit  measurement  changes  are 
also  reflected  in  the  diagram  of  Figure  4. 

To  complete  the  requirements  of  the  EASY  modelling 
program,  a component  block  diagram  with  interconnections  was 
specified.  The  output  of  the  EASY  Modelling  Generation 
program  provided  a schematic  of  the  overall  system  indicating 
input,  output,  and  parameter  requirements  of  the  flight 
control  system  aircraft  model  and  pilot  model.  A listing 
of  the  model  description  statements  along  with  a computer 
generated  schematic  diagram  of  the  total  system  modelled 
is  included  in  Appendix  A. 

F-16  EASY  Analysis 

Parameter  requirements  of  the  analysis  program  for  the 
flight  control  system  and  pilot  model  were  specified  by  the 
input  list  generated  by  the  model  program.  This  data  is 
in  two  parts;  first,  data  necessary  to  generate  the  table 
look-up  functions  of  the  flight  control  system;  and  secondly, 
parameter  values  of  the  standard  components.  The  gain 
scheduling  blocks  indicated  in  Figure  2 were  built  using 
the  table  function  components,  FU.  This  tabular  data  was 
loaded  by  describing  both  the  independent  and  dependent 
variables.  Additional  parameters  were  loaded  in  the  analysis 
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program  to  specify  either  linear  interpolation  or  extrapola- 
tion of  the  table  look-up  functions.  Following  the  tabular 
data,  parameter  values  were  loaded  to  specify  all  require- 
ments of  the  standard  components  shown  in  Figure  3 and 
Figure  4. 

Additional  data  requirements  included  stability  deriva- 
tive information  necessary  to  satisfy  the  longitudinal  equa- 
tions of  motion.  Inconsistencies  were  found  in  stability 
derivative  definitions,  sign  conventions,  and  unit  measure- 
ments. It  was  therefore  necessary  to  develop  an  axis  system 
and  sign  convention  to  be  consistent  throughout  the  simula- 
tions. The  most  convenient  system  was  a set  of  mutually 
perpendicular  reference  axes  intersecting  at  the  center  of 
gravity  of  the  aircraft.  About  this  point,  the  aircraft 
motion,  moments  and  forces  were  measured.  The  positive 
directions  for  these  axes  were  selected  as:  forward  or 
opposite  the  direction  of  airflow  for  the  X axis;  to  the 
right  for  the  Y axis;  and  downward  for  the  Z axis.  Reference 
data  provided  by  the  aircraft  manufacturer  and  the  test  case 
from  the  Griffin  program  was  selected  using  the  stability 
axes  as  reference.  The  stability  axes  are  established  with 
the  X axis  parallel  to  the  undisturbed  airflow  with  respect 
to  the  aircraft  body.  The  axis  system  selected  is  described 
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in  Figure  5.  The  sign  convention  established,  although 

not  universally  used,  is  consistent  with  both  the  data 

presented  by  General  Dynamics  and  that  of  the  test  case 

selected.  A graphic  description  of  the  sign  convention 

used  is  shown  in  Figure  6 (Ref  11). 

Data  preparation  for  the  constant  coefficient  aero 

model  of  the  EASY  analysis  program  included  external  forces, 

torques,  and  aerodynamic  stability  derivatives.  The  non- 

■ 

dimensional  stability  derivatives  were  obtained  from  test 
data  derived  from  a test  program  used  by  the  Air  Force 
Flight  Dynamics  Laboratory  LAMARS  simulation  program.  Run 
| #43  of  this  Griffin  program  was  used  as  a data  base  and  is 

listed  in  Appendix  B.  It  should  be  noted  that  the  dynamic 
stability  derivatives  (e.g.  functions  of  control  surface 
deflections)  are  listed  in  per  degree  units  while  static 
stability  derivatives  are  listed  in  radian  measure.  The 
non-dimensional  derivatives  used  to  satisfy  the  longitudinal 
equations  of  motion  of  the  EASY  program  are  shown  in  Table 
III  along  with  their  respective  EASY  program  names  (Ref  6). 

F-16  Aircraft  Model  Validation 

With  the  tabular  data,  parameter  values,  and  stability 
derivative  requirements  satisfied,  the  EASY  analysis  program 
was  used  to  verify  the  aircraft  dynamics.  The  first  check 
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Figure  5.  Axes  System  Diagram  (Ref  11) 


Figure  6.  Simulation  Sign  Convention  Diagram  (Ref  11) 
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of  the  system  was  to  trim  the  aircraft  in  straight  and  level 
flight  at  an  altitude  of  20,000  feet  and  airspeed  of  .8 
Mach,  or  an  equivalent  829.5  feet  per  second.  In  order  to 
trim  the  aircraft,  two  additional  integrators  were  added 
to  the  system  model.  The  input  of  the  first  trim  integrator 
was  the  difference  in  reference  altitude  of  20,000  feet  and 
the  actual  calculated  altitude  of  the  SD  component  during 
the  trim  iterations.  This  difference  error  was  then  inte- 
grated and  the  output  fed  into  the  system  as  a stick  input 
to  help  achieve  a steady  state  at  20,000  feet.  This  trim 
integrator  is  shown  in  the  upper  left  hand  corner  of  Figure 
3.  Likewise,  a second  trim  integrator  was  used  to  compare 
actual  velocity  calculated  in  the  SD  component  to  the 
reference  velocity  of  829.5  feet  per  second.  This  velocity 
difference  was  integrated  and  the  error  was  used  to  help 
achieve  the  desired  trim  condition. 

Through  the  use  of  program  commands,  a steady  state 
system  solution  was  determined.  Following  trim  iterations 
of  the  EASY  program,  the  computer  output  shown  in  Figures 
7 and  8 verify  the  aircraft  trim  condition.  The  system 
states,  as  defined  by  the  EASY  program,  are  those  quantities 
described  by  first  order  differential  equations.  As  shown 
in  Figure  7,  the  system  modelled  consists  of  25  states.  The 


Figure  7.  EASY  Program  Steady  State  Output 
(States,  Rates,&  Variables) 
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Figure  8.  EASY  Program  Steady  State  Output 
( Parameters ) 


value  of  each  stale  is  given  for  the  trim  condition.  Out- 
put variables  are  shown  that  correspond  t:o  the  component 
outputs  of  Figure  3,  p.  21.  The  parameter  values  for  each 
standard  component  are  listed  in  Figure  8.  To  achieve  the 
desired  trim  condition,  the  pilot  model  was  isolated  from  the 
aircraft  and  flight  control  models.  Program  commands  were 
then  used  to  generate  steady  state  iterations.  In  doing 
this,  lire  dynamic  equations  of  motion  were  perturbed  after 
each  iteration  step  until  the  trim  condition  specified  by 
altitude  and  airspeed  initial  conditions  was  achieved.  The 
state  variables  indicate  the  quantities  that  resul t from 
integrating  the  set  of  first  order  differential  equations 
that  comprise  the  dynamic  system  model. 

Once  that  an  operating  point  was  established,  all  initial 
conditions  of  integrator  states  were  transferred  to  the  sys- 
tem through  a program  command.  Witli  the  two  additional  trim 
integrators  for  altitude  and  airspeed  turned  off  (i.e.  inte- 
grator states  "frozen")  verification  of  the  F-16  aircraft 
characteristics  continued. 

The  model  developed  by  the  KASY  program  which  includes 
the  K-16  aircraft  character is  tics,  longitudinal  flight  con- 
trol system,  and  longitudinal  pilot  model,  is  described  by 
the  simplified  block  diagram  of  Figure  9. 
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Figure  9.  EASY  Model  Simplified  Diagram 

By  selecting  summing  junction  parameters  to  open  all 
flight  control  system  feedback  loops  and  freezing  integra- 
tors in  the  pilot  model  and  flight  control  system,  transfer 
functions  of  the  aerodynamic  variables  for  elevator  deflec- 
tion could  be  obtained.  With  the  above  provisions  completed 
through  program  commands,  the  total  system  was  effectively 
reduced  to  the  aircraft  dynamic  model  shown  in  the  simpli- 
fied diagram  of  Figure  10.  The  basic  aircraft  dynamics  of 
the  Griffin  program  were  compared  with  EASY  Analysis  program 
calculations.  A comparison  of  the  computer  generated 
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Figure  10.  EASY  Transfer  Function  Diagram 


character istic  equation  roots  with  the  Griffin  program  data 
is  listed  in  Table  IV. 


Table  IV 

Roots  of  the  F-16  Characteristic  Equation 
(Longitudinal  Dynamics) 


Griffin  Data 

Easy  Analysis  Data 

_____ 

Short  Period 

+.8282 

+.8168 

-2.745 

-2.746 

Phugoid 

-.0262  + j - 1528 

-.0288  + j .1245 

-.0262  - j . 1528 

-.0288  - j . 1245 

At  the  flight  condition  selected  for  this  simulation, 
flight  data  indicates  that  the  basic  F-16  airframe  is  un- 
stable. The  relaxed  static  stability  is  demonstrated  by 
the  real  positive  eigenvalue.  As  shown  above,  the  short 
period  mode  of  the  F-16  aircraft  has  two  real  roots.  One 
of  these  in  the  right  half  plane  brings  about  the  airframe 
instability.  The  phugoid  mode  is  indicated  by  the  pair  of 
dominant  complex  poles  near  the  origin. 

Although  the  numerator  characteristics  were  calculated 
in  the  EASY  Analysis  program,  it  was  not  possible  to  have 
these  printed  to  enable  a comparison  with  the  numerator 
dynamics  of  the  Griffin  program.  However,  transfer  functions 
were  calculated  using  the  EASY  program  commands  and  the 
computer  generated  Bode  diagrams  were  used  to  validate  the 
numerator  dynamics.  From  the  Griffin  program  data,  the 
following  analytical  transfer  functions  of  angle  of  attack 
per  elevator  deflection  and  pitch  rate  per  elevator  deflec- 
tion were  derived: 

= - . 1329 ( s + 148.5 ) ( s^  + ,0177s  + .0040)  (3) 

(s  - .8282) ( s + 2 .745 ) ( s2  + .0524s  + .0240) 

= -19.72s(s  + . 01741 )( s + 1.312) (4) 

(s  - . 8282 ) ( s + 2.745 )( s2  + .0524s  + .0240) 
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The  AFIT  frequency  response  program  (FREQR)  was  used 
to  generate  plots  of  the  above  transfer  functions.  Figure 

11  shows  the  angle  of  attack  per  elevator  deflection  trans- 
fer function  of  the  Griffin  program  data.  This  compares 
very  favorably  with  the  computer  generated  EASY  magnitude 
and  phase  plots  of  alpha  per  delta  elevator  shown  in  Figures 

12  and  13,  respectively.  A similar  comparison  of  the  pitch 
rate  per  delta  elevator  transfer  functions  was  made  and  the 
results  are  shown  in  Figures  14-16.  The  above  comparisons 
clearly  serve  to  substantiate  the  validation  of  the  F-16 
aircraft  model  developed. 

As  mentioned  earlier,  the  F-16  basic  aircraft  is  un- 
stable at  the  flight  condition  selected  for  this  simulation. 
This  instability  is  evidence  of  a positive  static  margin 

which  is  defined  as  the  ratio  of  C to  C.  . The  instability 

m L J 

a a 

of  the  basic  F-16  airframe  was  demonstrated  by  using  the 
simulation  program  commands  of  the  EASY  Analysis  program. 

With  all  integrator  states  in  the  flight  control  system 
frozen  and  also  insuring  that  all  feedback  channels  of  the 
flight  control  system  were  opened,  a time  simulation  was  run 
to  show  the  dynamics  of  the  aircraft  alone  with  no  flight 
control  system.  An  initial  condition  equivalent  to  a 
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FREQUENCY,  RPS 

Figure  11.  Alpha/Delta  Elevator  Transfer  Function  of  the  Griffin  Program  Data 


42 


Figure  12.  EASY  Generated  Magnitude  Plot  of  Alpha/Delta  Elevator  Transfer  Function 


FREQJENCY,  RPS 

Rate/Delta  Elevator  Transfer  Function  of  the  Griffin  Program  Data 


1 

one  degree  angle  of  attack  pertubation  was  input  into  the 
aircraft  system  and  the  resulting  dynamic  response  is  shown 
in  Figure  17.  It  can  be  seen  that  the  aircraft  does  not 
return  to  its  equilibrium  trim  condition  with  even  a slight 
angle  of  attack  pertubation.  The  display  shows  the  angle  of 
attack,  pitch  angle,  and  altitude  increasing  while  the  air- 
craft airspeed  decreases.  This  demonstrates  the  necessity 
of  maintaining  the  flight  control  system  to  harness  the 
inherent  aerodynamic  instability. 

EASY  System  Analysis 

The  EASY  Analysis  program  was  further  used  to  examine 
the  present  F-16  flight  control  system.  Investigations 
continued  to  establish  a measure  of  the  system  effectiveness. 

The  present  F-16  flight  control  system,  which  is  predomin- 
antly normal  acceleration  feedback,  was  evaluated  in  terms 
of  tracking  performance.  As  shown  in  Figure  18,  a closed 
loop  system  was  established  by  having  the  pilot  model  respond 
to  an  angular  error  input.  The  command  angle  was  the  dif- 
ference between  the  aircraft  pitch  angle  and  a prescribed 
reference  angle.  By  establishing  this  pseudo  tracking  task, 
a first  order  approximation  of  the  director  sight  implementa- 
tion was  achieved.  The  pilot  model  parameters  were  selected 
to  be  consistent  with  the  director  sight  characteristics. 
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Figure  18.  EASY  Closed  Loop  System  Schematic  Diagram 


Control  of  the  pitch  feedback  summing  junction  of  Figure  18 
was  gained  through  EASY  program  commands  and  by  selecting 
the  appropriate  parameter,  the  0/0p£ p transfer  function  was 
established  for  either  closed  or  open  loop  analysis. 

The  magnitude  and  phase  Bode  plots  of  the  normal  accel 
eration  configuration  for  the  closed  loop  are  shown  in 
Figures  19  and  20,  respectively.  The  magnitude  plot  indi- 
cates the  transfer  function  is  well  behaved  and  the  maximum 
peak,  M^,  of  1.43  occurs  at  a frequency  of  .356  radians  per 
second.  As  seen  in  Figure  19,  the  system  bandwidth,  is 


approximately  .7  radians  per  second.  Next  the  open  loop 
transfer  function  was  investigated  by  opening  the  feedback 
loop  at  the  summing  junction.  The  Bode  plots  of  the  open 
loop  normal  acceleration  transfer  function  are  shown  in 
Figures  21  and  22.  As  seen  in  Figure  21,  the  gain  cross- 
over occurs  at  approximately  .5  radians  per  second  and  the 
phase  margin  (Figure  22)  is  52  degrees.  To  evaluate  the 
effectiveness  of  the  pseudo  tracking  task  of  the  pilot  model, 
a time  simulation  was  run  to  show  the  system  response  to  a 
step  input  of  5.5  degrees  pitch.  The  time  response  to  this 
step  input  is  shown  in  Figure  23. 

The  EASY  Analysis  program  was  also  used  to  investigate 
the  merits  of  a pitch  rate  command  system.  As  discussed  in 
Chapter  I,  the  C*  design  concept  fortheF-16  flight  control 
system  employs  a predominantly  normal  acceleration  command 
system  at  cruise  airspeeds.  This  is  evidenced  by  the  larger 
weighting  of  normal  acceleration  to  pitch  rate  in  the  feed- 
back channels  of  the  longitudinal  control  system.  For  air 
combat  tracking  tasks,  pitch  rate  feedback  could  be  made 
predominant  by  adjusting  this  weighting  relationship.  One 
possible  implementation  would  be  the  total  weighting  of  pitch 
rate  as  the  command  signal  with  the  elimination  of  normal 
acceleration  feedback. 
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FREQUENCY,  RPS 

Figure  22.  Normal  Acceleration  Configuration  Open  Loop  Bode  Phase  Plot 
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Figure  23.  Normal  Acceleration  Configuration  Time  Response  to  Step  Input  Tracking  Error 


This  method  was  attempted  and  the  EASY  program  was 
modified  to  incorporate  a pitch  rate  control  system.  As 
shown  in  the  model  diagram  of  Figure  3,  p.  21,  the  normal 
acceleration  channel  was  eliminated  as  a feedback  signal  by 
setting  the  gain  parameter  equal  to  zero.  The  angle  of 
attack  feedback  channel  was  maintained  to  aid  stability. 
Because  the  primary  feedback  variable  was  to  be  pitch  rate, 
the  wasnout  filter  in  the  pitch  rate  channel  was  removed. 

The  lead- lag  transfer  function  2(s  + 15 ) / ( s + 30)  was  im- 
plemented as  the  new  LEE3  component  to  increase  the  system 
bandwidth.  Additionally,  a root  locus  analysis  was  used  to 
determine  that  a pitch  rate  feedback  system  gain  of  .3  would 
provide  an  overall  system  damping  factor  of  approximately  .7. 
An  improved  system  performance  was  indicated  by  evaluating  the 
open  and  closed  loop  transfer  functions.  The  system  response 
with  the  pitch  rate  feedback  system  is  indicated  in  Figures 
24,  25,  26,  and  27.  The  response  of  the  pitch  rate  system  is 
similar  to  that  of  the  normal  acceleration,  however,  certain 
points  are  noteworthy.  For  example,  the  maximum  peak  of  the 
closed  loop  transfer  function,  shown  in  Figure  24,  occurs  at 
the  same  frequency  as  in  the  normal  acceleration  system,  but 
reduced  to  1.28.  As  noted  also  in  Figure  24,  the  effective 
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bandwidth  of  the  system  has  been  increased  from  .7  to  1.0 
radians  per  second.  Figure  26  indicates  the  gain  cross- 
over of  the  open  loop  pitch  rate  transfer  function  occurs 
at  approximately  .75  radians  per  second. 

The  above  observations  imply  that  a pitch  rate  command 
system  could  be  successfully  implemented  to  achieve  overall 
system  performance  improvement.  To  verify  this,  the  time 
response  of  the  pitch  rate  system  to  a step  input  tracking 
error  is  shown  in  Figure  28.  The  faster  response  with  less 
overshoot  is  evidence  of  an  improved  system. 

To  summarize  the  results  of  the  normal  acceleration  and 
pitch  rate  control  system  investigation.  Table  V lists  a 
comparis  m of  quantities  of  interest  (Ref  14). 

Pilot  Model  Study 

The  pilot  model  adapted  for  this  investigation  was 
developed  by  McDonnell  Douglas  Corporation  with  gain  para- 
meters adjusted  by  General  Dynamics  to  meet  the  character- 
istics of  the  F-16  aircraft.  It  was  beyond  the  scope  of 
this  thesis  to  extensively  evaluate  the  pilot  model  and  con- 
duct an  in-depth  study  to  confirm  that  the  pilot  model  used 
is  adequate  for  the  F-16  aircraft  and  gunsight  characteris- 
tics. However,  a brief  examination  of  the  pilot  model  loop 
was  done  to  insure  stability  and  adequate  performance  with 
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Table  V 

Quantities  of  Interest 
Normal  Acceleration  vs  Pitch  Rate  Feedback 

(Ref  14) 


Normal  Acceleration 
Feedback 

Pitch  Rate 
Feedback 

Max  Peak  Value 

M 

m 

1.43 

1.28 

Peak  Frequency 

id  (rps) 
m 

.356 

.356 

F.ase  Margin 
<j>c  ( deg ) 

52 

63  ' 

Gain  Margin  Frequency 

a)  ( rps ) 
c 

.52 

.75 

Gain  Margin 
(dB) 

22.2 

18.4 

Bandwidth 

(rps) 

.7 

1.0 

Time  Simulation  Results 


Peak  Overshoot 


M0  (%) 

17.3 

13.6 

Peak  Time 

T (sec) 
p 

4.7 

3.25 

Settling  Time 

T (sec) 
s 

9.5 

7.2 

Rise  Time 

Tr( sec) 

.5 

.4 

Min 


|w  jn 


* 


i 

j 
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the  pitch  lute  control  system.  For  the  longitudinal 
the  pilot  loop  included  Hie  long  i 1 udinal  pilot  model 


shown  in  Fipui  e 1,  p.  10,  as  well  a;-,  the  pitch  command  stick 
gradient  and  the  stick  conditioning  lag  filter.  These  three 
components  represent  the  pilot  loop  and  the  open  loop  trans- 
fer function  is  as  follows: 


C(s) 


119.1 (s3  r 1 ,32bs"  + 1 .lbs  r 1 ■ 1 25 ) 
s ( s + 8 . 8 ) ( s + 20 . ) ( s'  + 1 . 2s  + 1 . ) 


(5) 


where  the  pitch  command  stick  gradient  gain  was  selected 
as  the  upper  slope  value  of  .8182  as  shown  in  the  Holdout 
diagram  of  Figure  2,  p.  20. 

Again,  the  AFIT  FREQR  program  was  used  to  produce  the 
Bode  magnitude  and  phase  plot  of  Figure  29.  The  quadratic 
low  pass  filter  of  the  pilot  model  causes  an  extensive  phase 
shift  near  a frequency  of  1.  radian  per  second.  An  accept- 
able phase  margin  .is  maintained  up  to  frequencies  near  10 
radians  per  secoiTd.  The  magnitude  plot,  although  well 
behaved  throughout  the  operational  region  of  the  pilot  model, 
i ndic.it . • de.drabl.  - 20  dB  pei  decade  slope  at  the  lower 

i tin  •’  s:  \ <•!  ! vpe  model  representation  of 

•r  * i ' . ITi  i l : : ‘ i . i t o:  that  no  adverse 
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effects  should  be  expected  with  the  same  pilot  model  im- 
plemented in  the  pitch  rate  system  since  the  effective  band- 
width of  the  system  has  been  increased  from  .7  to  1.0 
radian  per  second.  In  fact,  an  improvement  may  be  attained 
since  larger  magnitude  outputs  may  be  achieved  through 
operation  at  the  slightly  higher  frequencies  possible  with 


the  pitch  rate  control  system. 


EASY  Analysis  Summary 


Verification  of  the  F-16  aircraft  dynamics  indicated 
that  the  EASY  model  characteristics  were  consistent  with  the 
simulation  test  data.  By  implementing  an  analytical  pilot 


model,  closed  loop  simulation  analysis  was  made  possible. 

A pseudo  target  tracking  system  was  developed,  and  the 
performance  of  the  present  F-16  flight  control  configuration 
was  examined.  Investigation  of  a proposed  pitch  rate  flight 
control  configuration  provided  both  frequency  domain  and 
time  domain  results  that  indicated  target  tracking  improve- 
ments were  possible  by  implementing  a pitch  rate  control 
scheme.  A listing  of  the  EASY  program  statements  used  to 


generate  the  control  system  analysis  is  given  in  Appendix  D. 


V.  Development  of  a Simulation  Program 


TAWDS  Program  Discussion 

To  evaluate  the  air-to-air  tracking  performance  of  the 
F-16  aircraft  model,  a simulation  program  was  needed.  The 
Terminal  Aerial  Weapon  Delivery  Simulation  (TAWDS)  program 
produced  by  McDonnell  Douglas  Corporation  (Ref  15)  was  well 
suited  to  simulate  the  air-to-air  encounters  and  provide  an 
evaluation  of  the  flight  control  systems  employed.  Although 
the  program  has  provisions  to  simulate  weapon  delivery  tasks 
for  air-to-air  gunnery,  air-to-ground  gunnery,  and  bombing, 
only  the  air-to-air,  aerial  gunnery  programs  were  included 
in  this  study.  The  TAWDS  digital  simulation  program  was 
used  to  simulate  both  the  present  F-1S  flight  control  con- 
figuration and  the  pitch  rate  configuration  of  Chapter  IV. 
Implementing  the  pilot  model  discussed  in  Chapter  II,  a 
deterministic  evaluation  of  the  aircraft's  tracking  perform- 
ance was  completed.  Provisions  of  the  TAWDS  program  include 
many  factors  associated  with  aerial  gunnery.  For  example, 
the  program  includes  the  modelling  effects  of  aircraft 
dynamics,  control  system  characteristics,  gunsight  character- 
istics, pilot  control  parameters,  attacker  to  target  geometry, 
target  maneuvering,  gun  orientation,  gun  rate  of  fire  and 
recoil  forces,  bullet  trajectories,  random  windgusts,  and 


stationary  source  errors.  The  above  considerations  make  the 


TAWDS  digital  simulation  program  a well  suited  analytical 


tool  to  evaluate  the  non-linear  six-degree-of -freedom  air- 


to-air  terminal  tracking  task. 


TAWDS  Programming  Techniques 


The  deterministic  mode  of  the  TAWDS  air-to-air  program 


uses  non-linear  time  varying  equations  to  simulate  a six- 


degree-of-  freedom  attacking  aircraft  tracking  and  firing  at 


a f ive-degree-of-f reedom  maneuvering  target.  The  major 


subroutines  of  the  TAWDS  air-to-air  program  are  called  by 


the  Executive  subroutine  to  describe  the  air-to-air  terminal 


weapon  delivery  task.  These  subroutines  include  Data  Input, 


Initial  Encounter,  Initial  Condition,  Measurement  Error, 


Airframe,  Augmentation,  Pilot,  Target  Initialization,  Target 


Aircraft,  Relative  Geometry,  Bullet  Time  of  Flight,  LCOS 


Sight,  Director  Sight,  Bullet  Integration,  Performance, 


Runge-Kutta  Integration,  and  Output  Subroutines  (Ref  15). 


In  preparation  for  using  the  TAWDS  program,  it  was 


necessary  to  develop  tabular  data  for  the  six-degree-of- 


freedom  non-linear  F-16  aircraft  model.  Extensive  stability 


derivative  data  supplied  by  the  aircraft  manufacturer  (Ref 


16)  was  input  into  the  program  to  provide  table  look-up 


parameters  necessary  to  satisfy  the  six-degree-of-f reedom 
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equations  of  motion.  A flight  condition  of  altitudes  near 
20,000  feet  with  airspeed  varying  from  .8  to  .9  Mach  was 
considered.  The  aircraft  model  selected  as  the  data  base 
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for  the  EASY  programming  was  again  used  in  the  TAWDS  pro- 
gram. Therefore,  the  aircraft  characteristics  of  Table  III, 
p.  31,  were  implemented.  Additional  simulation  specifica- 
tions were  allowed  and  the  reader  is  referred  to  Appendix  C 
and  Ref  15  if  programming  details  are  desired. 

The  generic  design  of  the  TAWDS  Data  Input  subroutine 
allowed  easy  implementation  of  parameter  values.  Different 
aircraft  characteristics,  flight  control  or  pilot  model 
parameters,  or  gunsight  selections  could  be  made  through 
data  changes.  Initially,  use  of  the  TAWDS  program  was  to  be 
limited  to  the  longitudinal  axis  of  the  flight  control  sys- 
tem, however,  since  weapon  system  effectiveness  was  the 
overall  objective  of  the  study,  it  did  not  seem  realistic 
to  evaluate  this  system  in  only  one  plane  of  motion.  It 
was  hoped  that  air-to-air  combat  encounters  could  be 
developed  to  provide  a tracking  task  that  would  realistically 
evaluate  the  candidate  flight  control  systems.  Therefore, 
it  was  decided  to  implement  a full  six-degree-of-freedom 
simulation. 
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After  the  aircraft  data  requirements  were  satisfied, 
the  longitudinal  flight  control  system  was  implemented. 

Since  the  F-16  has  a very  non-conventional  flight  control 
system,  the  generic  flight  control  models  of  the  TAWDS  pro- 
gram were  not  adaptable  for  the  F-16  aircraft.  Therefore, 

FORTRAN  statements  were  used  instead  to  develop  flight 
control  signal  processing.  The  same  control  system  simpli- 
fications in  the  longitudinal  axis  as  employed  in  the  EASY 
programming  were  also  considered  in  the  TAWDS  program.  Since 
the  primary  objective  was  to  evaluate  longitudinal  flight 
control  tracking  characteristics,  the  cross  coupling  of  lat- 
eral and  longitudinal  control  signalling  was  also  eliminated. 

FORTRAN  statements  were  used  to  generate  the  gain  scheduling 
requirements  of  the  flight  control  system.  Functions  such 
as  dynamic  pressure  adjusted  for  compressibility,  static 
pxessure,  and  Mach  number  variables  were  developed  within 
the  logic  of  the  flight  control  subroutine.  The  schematic 
diagram  shown  in  Figure  30  shows  the  longitudinal  flight 
control  system  implemented  in  the  TAWDS  program.  The  lat- 
eral-directional flight  control  system  for  the  F-16  model  is 
shown  in  Figure  31.  The  parameter  names  and  transfer  function 
names  as  indicated  in  the  block  diagram  hold  no  significance 
except  they  were  programmed  parameters  and  transfer  function 
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Figure  31.  Continued 


names  already  used  in  the  generic  TAWDS  program.  These 
were  again  used  in  the  FORTRAN  programming.  Using  these 
pre-programmed  names  simplified  programming  efforts,  but 
these  names  bear  no  relation  to  the  generic  TAWDS  program. 

An  additional  trim  subroutine  was  also  vsed  to  replace 
the  aircraft  trimming  techniques  of  the  TAWDS  program.  This 
subroutine  perturbed  the  six-degree-of-f reedom  equations  by 
varying  aircraft  angle  of  attack,  forces,  and  moments  until 
a satisfactory  trim  solution  was  achieved  about  the  desired 
flight  condition.  An  added  feature  allows  specification  of 
the  trim  condition  in  terms  of  g's  by  reading  in  the  para- 
meter name  GKMECH.  For  example,  to  achieve  a trim  condition 
of  straight  and  level  unaccelerated  flight,  GKMECH  = 1.0. 
After  selecting  a 1 g flight  condition,  the  TAWDS  program 
output  of  stability  derivatives  and  trim  values  indicated 
only  minor  numerical  variance  whs  it  compared  ho  hhs  hr  mi 
condition  values  of  the  EASY  program. 

Next  it  was  necessary  to  initialise  the  flight  control 
system.  The  values  of  the  aircraft  feedback  variables  were 
used  to  initialise  the  transfer  functions  of  the  flight  con- 
trol system.  The  first  pass  through  the  flight  control  sub- 
program (AUTSI1)  initialised  the  transfer  functions  and  set 
control  conditions  to  satisfy  the  trim  condition.  A stick 
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trim  bias  as  shown  in  the  pilot  channel  of  Figure  30  was 
calculated  to  provide  the  nominal  stick  forces  necessary  to 
balance  the  control  surface  deflection. 

The  pilot  model  and  gain  parameters  as  discussed  in 
Chapter  III  were  also  included  in  the  TAWDS  programming 
model.  Since  a full  six-degree-of-freedom  simulation  was  to 
be  developed,  the  multi-axis  F-16  pilot  model  was  adapted 
which  included  not  only  the  longitudinal  axis,  but  also 
the  lateral-directional  axis. 

The  generic  quality  of  the  TAWDS  programming  ability 
was  very  helpful  in  implementing  the  pitch  rate  control  sys- 
tem of  Chapter  IV.  It  was  not  necessary  to  extensively 
change  the  flight  control  model  as  in  the  EASY  program. 
Parameter  value  changes  read  in  as  data  could  effectively 
modify  the  normal  acceleration  system  to  incorporate  the 
pitch  rate  control.  For  example,  setting  the  parameter 
GKNZ2  shown  in  Figure  30  equal  to  aero  eliminated  normal 
acceleration  feedback.  The  transfer  function  CBQQ/CCQQ  was 
set  to  unity  and  the  transfer  function  CPCOML/CCCOML  was 
changed  to  the  desired  filter  2(s  + 15)/(s  + 30).  By  set- 
ting the  parameter  GKF  = .3  as  determined  by  the  root  locus 
analysis  mentioned  in  Chapter  IV,  and  adjusting  the  pilot 
loop  gain  to  command  degrees  per  second  instead  of  g's, 
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(e.g.  CBSE  = 10.00),  the  pitch  rate  system  implementation 
for  the  TAWDS  program  was  completed. 


Data  parameters  would  also  allow  the  selection  of  either 
the  LCOS  or  director  sight  to  be  implemented  into  the  simula- 
tion program  for  each  run.  By  implementing  either  of  the 
above  sights,  the  system  developed  to  this  point  could  be 
considered  as  a manual  director  system  or  a manual  LCOS 
system.  Variations  of  the  multi-axis  pilot  parameters  would 
be  necessary  to  operate  the  different  gunsight  systems. 

Aii-to-Air  Scenarios 

A series  of  air-to-air  gunnery  encounters  were  developed 
to  evaluate  the  flight  control  systems.  The  Initial  Condi- 
tion subroutine  was  programmed  to  set  the  attacker  aircraft 
at  the  flight  condition  to  begin  eight  air-to-air  engagements. 
For  each  target  scenario,  the  attacker  aircraft  was  initia- 
lized near  straight  and  level  unaccelerated  flight  at 

20.000  feet  and  airspeed  of  .8  Mach.  A target  range  of 

5.000  feet  was  selected  for  the  first  four  target  maneuvers 
and  a range  of  7,000  feet  for  the  last  four  maneuvers.  Vary- 
ing range  rates  (i.e.  attacker  to  target  closure  rate)  from 
-50  to  -200  feet  per  second  were  programmed  also.  This 
established  the  initial  target  aircraft  airspeed  for  each 
scenario . 
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In  addition  to  initializing  the  attacker  aircraft,  tar- 
get aircraft  maneuvers  were  developed.  To  be  consistent 
with  terminal  tracking  environments,  evasive  maneuvers 
were  programmed  for  an  F-4  aircraft  target  model.  Speci- 
fications for  the  target  model  included  time  intervals  to 
begin  and  end  each  maneuver,  bank  angle  rate,  angle-of- 
attack  rate,  and  thrust  rate.  A three  digit  code  of  a 
selected  bank  angle  rate,  angle-of-attack  rate,  and  thrust 
rate  was  developed  for  specific  target  maneuvers  at  selected 
time  intervals  during  the  tracking  scenario.  Table  VI 
describes  the  maneuvers  developed  and  the  corresponding  NIC 
number  used  to  select  each  encounter.  Generally,  the  tar- 
get maneuvered  during  the  first  two  seconds  of  the  program 
into  a right  or  left  90°  bank  turn  for  a specified  g and 
then  completed  either  a split  S or  jink  maneuver  at  a 
specified  time  as  the  scenario  reached  the  terminal  time 
of  fifteen  seconds.  These  eight  air-to-air  encounter 
scenarios  were  used  as  a basis  to  evaluate  the  director 
sight  implementation  under  both  the  normal  acceleration 
and  pitch  rate  flight  control  configurations.  With  the 
deterministic  mode  selected,  the  program  would  complete  one 
pass  through  each  case  and  generate  time  history  data  to 
measure  the  attacker  aircraft  performance  while  tracking 
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the  target  aircraft  through  the  selected  maneuvers.  Both 
mean  and  RM3  elevation  and  traverse  tracking  error  informa 
tion  was  to  be  used  to  measure  the  tracking  effectiveness. 
The  TAWDS  program  information  is  shown  in  Appendix  C. 


Table  VI 

TAWDS  Air-to-Air  Target 


Maneuvers 


— ■ — 1 

NIC 

Max 

g 

Closure 
( f t /sec ) 

Time 
( sec ) 

Range 

(ft) 

Maneuver 

Description 

1 

3.6 

-100 

0-2 

5000 

Roll 

90° 

left 

15 

1000 

2 

4.2 

-50 

0-2 

5000 

Roll 

90° 

left 

10-11 

2600 

Roll 

30° 

left 

15 

1200 

3 

3.2 

-150 

0-2 

5000 

Roll 

90° 

left 

12-13 

1500 

Roll 

30° 

right 

15 

800 

4 

3.4 

-150 

0-2 

5000 

Roll 

90° 

right 

. 

8-10 

2500 

Roll 

45° 

right 

15 

800 

5 

3.2 

-150 

0-2 

7000 

Roll 

90° 

right 

8-10 

4000 

Roll 

30° 

right 

15 

1400 

6 

3.5 

-200 

0-2 

7000 

Roll 

90° 

left 

8-10 

4000 

Roll 

45° 

left 

15 

1400 

7 

3.3 

-150 

0-2 

7000 

Roll 

90° 

left 

6-8 

5000 

Roll 

45° 

left 

15 

1400 

8 

4.3 

-200 

0-2 

7000 

Roll 

90° 

right 

6-8 

5000 

Roll 

90° 

right 

15 

1500 

VI.  Results 


TAWD 6 Simulation  Techniques 

The  results  of  the  EASY  Analysis  of  Chuptei  IV  ii.uica- 
ted  an  improved  target  tracking  system  could  be  achieved 
by  implementing  a pitch  rate  flight  control  system  for  the 
F-16  aircraft.  To  verify  these  results,  each  of  the  eight 
air-to-air  target  encounters  was  run  using  the  TAWD6  pro- 
gram model.  This  digital  simulation  included  the  full  six- 
degree-  of-freedom  aircraft  model,  multi-axis  pilot  model, 
and  lateral-directional  as  well  as  longitudinal  flight  con- 
trol system.  Both  the  present  F-16  flight  control  configu- 
ration and  the  proposed  pitch  rate  control  system  were 
evaluated  using  the  director  sight  implementation. 

Runge-Kutta  integration  was  used  for  the  simulation 
using  a .05  second  iteration  step  size.  Although  a smaller 
step  size  would  provide  better  numerical  accuracy,  it  would 
also  require  greater  computer  time.  The  step  size  of  .05 
was  determined  to  be  the  maximum  size  for  the  simulation 
components  and  this  step  size  worked  well  for  the  normal 
acceleration  configuration.  However,  for  the  proposed 
pitch  rate  configuration,  this  step  size  was  not  sufficient. 
Modelling  of  a time  constant  of  1/30  seconds  could  not  be 
accomplished  with  a step  size  of  .05  seconds.  It  was 


necessary  to  reduce  the  step  increment  size  to  less  than 
the  largest  time  constant  of  the  system.  Test  runs  were 
made  using  step  sizes  of  .01  and  .001  seconds.  Very  little 
numerical  differences  were  noted  with  the  smaller  step  size 
and  therefore  the  significant  increase  in  computer  operation 
time  was  not  justified.  Since  the  normal  acceleration  con- 
figuration system  had  already  been  completed  with  a step 
size  of  .05  second,  it  was  desirable  to  use  the  same  step 
size  for  the  pitch  rate  configuration  also.  To  do  this,  it 
was  necessary  to  slightly  modify  the  lead-lag  compensator 
that  had  been  implemented  for  the  pitch  rate  system.  The 
transfer  function  CBC0ML/CCCOML  was  adjusted  from  a value 
of  2(s  + 15 ) / ( s + 30)  to  2(s  + 7.5)/(s  + 15).  This  allowed 
the  .05  step  size  to  be  maintained  and  provided  significant 
computer  cost  savings. 

TAWDS  Simulation  Results 

The  tracking  performance  of  the  F-16  attacker  aircraft 
for  each  of  the  eight  encounters  is  summarized  in  Table  VII. 
The  mean,  RM5,  and  standard  deviation  of  both  the  elevation 
and  traverse  aiming  error  is  shown  for  both  systems.  The 
bar  graphs  of  Figure  32  and  33  compare  the  RMS  elevation  and 
traverse  aiming  errors,  respectively.  Improvements  in  ele- 
vation tracking  can  be  seen  while  the  traverse  errors  remain 
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Table  VII 

Summary  of  TAWDS  Simulation  Tracking  Aim  Error 
with  Director  Sight  Implementation 
(measured  in  mils) 


Present  F-16 

Normal  Acceleration  Configuration 


Mean 

RMS 

Standard 

Deviat ion 

NIC 

Elevation 

Traverse 

Elevation 

Traverse 

Elevation 

Traverse 

1 

-22.2 

-17.8 

24.1 

19.8 

9.4 

8.7 

2 

-21.8 

9.3 

24.3 

14.0 

10.7 

10.4 

3 

-21.2 

7.7 

23.3 

13.7 

9.6 

11.4 

4 

-24.0 

-24.8 

26.2 

27.9 

10.4 

12.8 

5 

-20.8 

-16.7 

25.7 

25.1 

15.0 

18.8 

6 

-21.6 

14.9 

26.8 

20.2 

15.9 

13.6 

7 

-20.6 

11.1 

26.3 

20.1 

16.3 

16.8 

8 

-21.4 

-28.1 

27.5 

35.4 

17.3 

21.5 

lgur a l loi 

1 

-15.8 

-17.1 

18.8 

19.4 

10.2 

9.2 

2 

-15.8 

8.9 

19.1 

14.5 

10.8 

11.5 

3 

-15.3 

7.3 

17.7 

14.0 

9.1 

11.9 

4 

-19.0 

-23.9 

24.0 

27.2 

14.6 

12.9 

5 

-13.2 

-13.7 

17.8 

27.8 

12.0 

24.2 

6 

-13.2 

13.4 

17.1 

22.4 

10.8 

18.0 

7 

-12.9 

9.2 

17.1 

22.5 

11.2 

20.5 

8 

-13.0 

-25.7 

16.8 

36.3 

10.7 

25.5 

RMS  ELEVATION  AIMING  ERROR 


RMS  TRAVERSE  AIMING  ERROR 


PITCH  RATE  CONFIGURATION 


^ - -kTi'  i,  . - 


approximately  the  same.  The  results  indicate  that  the  ele- 
vation aiming  error  has  been  improved  by  an  average  of 
20%  for  the  first  four  encounters  that  were  initialised  at 
a range  of  bOOO  feet.  The  improvement  is  increased  to  an 
average  of  3b%  for  the  last  four  encounters  that  were 
initialised  at  a range  of  7000  feet.  The  least  improvement 
of  all  is  shown  in  encounter  4.  This  is  attributed  to  the 
target  maneuvering  at  the  most  critical  range  for  air-to- 
air  combat  of  2500  feet.  It  is  interesting  to  note  that  the 
tracking  performance  of  the  pitch  rate  system  shows  the 
greatest  improvement  in  encounter  8.  This  was  considered 
the  most  difficult  encounter  of  all  since  the  F-4  target 
aircraft  completed  a split-S  maneuver  during  the  lb  second 
terminal  chase.  This  result  indicates  that  the  pitch  rate 
system  of  the  longitudinal  axis  provides  a greater  improve- 
ment over  the  normal  acceleration  configuration  as  the 
difficulty  of  the  tracking  task  increases. 


VII.  Conclusions  and  Roconunendat ions 


Conclus ions 


From  the  previous  analysis  and  simulation  results,  it 
is  concluded  that  improvements  in  target  tracking  are 
possible  by  tailoring  the  flight  control  system  of  fighter 
aircraft.  Both  the  EASY  Analysis  and  the  TAWDS  simulation 
results  indicate  improvements  in  the  system  response  of  a 
pitch  rate  feedback  implementation  as  compared  to  the  normal 
acceleration  implementation  presently  employed  in  the  F-16 
aircraft . 

EASY  Analysis . The  EASY  Analysis  program  allowed  inves- 
tigations of  system  responses  with  both  the  normal  accelera- 
tion and  pitch  rate  flight  control  configurations.  By 
eliminating  the  normal  acceleration  feedback  and  compensa- 
ting the  pitch  rate  feedback,  it  was  shown  using  the  EASY 
program  that  a very  similar  overall  closed  loop  system  was 
maintained,  with  certain  noted  improvements.  A better 
damped  system  was  achieved  with  an  increased  phase  margin 
and  bandwidth  as  shown  by  Table  V,  p.  63.  The  time  simula- 
tion results  of  the  EASY  program  verified  the  improved 
system  response  to  a pseudo  tracking  task  in  the  longitu- 
dinal axis.  The  aircraft  mod*^  response  to  a pitch  angle 
step  input  to  the  pilot  model  indicated  improvements  with 
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the  pitch  rate  implementation.  As  shown  in  Table  V,  p.  63, 
the  pitch  rate  system  provides  better  damping,  faster  re- 
sponse, and  quicker  settling  to  the  step  input  signal  than 
the  normal  acceleration  system. 

TAW  PS  S imulat ion.  Similarly,  the  results  of  the  TAWDS 
simulation  air-to-air  combat  encounters  substantiated  the 
results  of  the  F.ASY  program  analysis.  The  tracking  aiming 
error  quantities  of  the  eight  air-to-air  encounters  as  shown 
in  Table  VII,  p.  S3,  clearly  indicate  the  improvements  of 
the  pitch  rate  flight  control  system.  Elevation  tracking 
error,  both  mean  and  RMS,  was  reduced  in  each  of  the  target 
tracking  scenarios.  These  air-to-air  encounters  provided 
a realistic  environment  to  test  the  performance  of  both  the 
normal  acceleration  and  pitch  rate  flight  control  configura- 
tions. The  different  target  airspeeds,  target  load  factors, 
and  target  rolling  maneuvers  of  each  scenario  allowed  the 
evaluation  of  each  system  throughout  a variety  of  flight 
conditions.  The  full  six-degree-of- freedom  non-linear 
simulation  developed  for  the  TAWDS  program  enabled  a realis- 
tic performance  evaluation  for  the  terminal  phase  of  air-to- 
air  combat. 

Aircra f t Model  Val idat ion.  Techniques  for  verifying 
the  aircraft  flight  dynamics  wore  shown  using  the  EASY 
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Analysis  program.  The  roots  of  the  <~haracter istic  equation 
for  the  longitudinal  dynamics  very  closely  matched  those  of 
the  LAMARS  test  case  used  to  model  the  F-16  aircraft  as 
shown  in  Table  IV,  p.  38.  In  addition,  the  transfer  func- 
tions of  alpha/elevator  deflection  and  theta  dot/elevator 
deflection  were  completed  to  verify  the  numerator  dynamics 
of  the  model . 

Pilot  Model  Performance.  The  satisfactory  performance 
of  the  analytical  F-16  pilot  model  clearly  demonstrated  the 
usefulness  of  such  a model  for  closed  loop  system  investi- 
gations. It  is  concluded  that  the  gain  sensitivity  para- 
meters of  the  pilot  model  are  sufficient  to  allow  a qualita- 
tive comparison  of  the  flight  control  configurations  ex- 
amined. The  pilot  model  control  stick  responses  of  the  TAWDS 
program  (see  Appendix  C)  indicated  reasonable  responses  to 
tracking  error  signals. 

Recommendat ions 

It  is  my  recommendation  that  extensive  digital  simula- 
tion be  continued  to  investigate  flight  and  fire  control 
systems  for  fighter  aircraft.  The  specific  recommendations 
from  this  study  include  the  following: 


Refine  Pitch  Rate  Control . The  mechanization  of  the 


flight  control  system  using  the  C*  concept  considered  only 
one  extreme  implementation,  that  of  eliminating  the  normal 
acceleration  and  conditioning  the  pitch  rate  feedback.  Since 
the  C*  approach  implies  a linear  blend  of  normal  accelera- 
tion, pitch  rate,  and  pitch  acceleration,  it  is  recommended 
that  other  possible  combinations  of  feedback  variables  be 
investigated. 

Flight  Envelope  Expansion.  The  end  game  problem  of 
air-to-air  combat  was  simulated  in  the  TAWDS  program.  The 
excursions  of  the  attacker  aircraft  did  not  effectively 
cover  the  operational  flight  envelope  of  the  F-16  aircraft 
during  the  15  second  encounters.  It  is  therefore  recommended 
that  more  extensive  flight  envelopes  be  considered. 

Pilot  Model  Validation.  The  analytical  pursuit  pilot 
model  used  in  this  thesis  was  based  on  actual  pilot  perform- 
ance data  during  target  tracking  tasks.  It  was  not  the  ob- 
jective of  this  study  to  verify  that  the  pilot  model  used 
closely  matched  that  of  the  actual  F-lb  pilot.  Instead,  it 
was  used  to  complete  a closed  loop  study  and  provide  con- 
sistent tracking  performance  for  the  different  flight  con- 
trol configurations  and  target  maneuvers.  If  it  is  required 
to  closely  match  the  actual  pilot  responses  for  future  F-16 
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simulations,  a detailed  study  of  the  pilot  model  is  recom- 
mended. Radical  departures  from  the  present  configuration 
would  indicate  the  need  for  a pilot  adaptation  study. 


Improve  Software  Programs . In  this  study,  the  EASY 
Model  Generation  and  Analysis  program  was  used  to  model  and 
analyze  the  present  flight  control  configuration  of  the  F-16 
aircraft  and  develop  a pitch  rate  implementation  to  improve 
target  tracking.  The  TAWDS  program  was  used  extensively  as 
a simulation  evaluation  tool.  However,  it  was  necessary  to 
master  the  programming  techniques  and  terminology  of  both 
programs.  Although  both  programs  provided  beneficial  re- 
sults, it  would  be  very  convenient  for  portions  of  each 
program  to  be  combined  in  one  compact  program.  The  EASY 
program  allowed  extensive  analysis  techniques  to  be  em- 
ployed by  simply  manipulating  parameter  values  and  state 
conditions  through  program  commands.  Both  frequency  domain 
and  time  domain  analysis  was  readily  available.  The  many 
subroutines  of  the  TAWDS  program  allowed  implementation  of 
many  of  the  considerations  for  air-to-air  combat.  The 
generic  quality  of  the  TAWDS  Input  subroutine  allowed  con- 
figuration changes  and  system  specifications  by  simply  man- 
ipulating data  values.  By  combining  the  benefits  of  each 
of  these  programs,  a very  large  reduction  in  the  time  and 
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effort  required  to  learn  both  programming  techniques  could 
be  achieved.  This  would  also  enable  the  user  to  develop  a 
greater  expertise  by  using  only  one  program.  In  particular, 
it  is  my  recominendat ion  that  the  frequency  response  and  time 
response  techniques  be  included  in  the  TAWDS  programming. 

In  this  manner,  a complete  program  package  could  be  main- 
tained so  that  the  programmer  could  readily  analyze  his 
implemented  system  and  complete  time  simulations  to  deter- 
mine system  effectiveness.  A reduction  in  computer  memory 
presently  required  to  run  both  the  EASY  and  TAWDS  programs 
would  significantly  reduce  the  operating  costs  and  improve 
program  turnaround  time. 

Apply  Techniques . The  programming  techniques  of  this 
thesis  in  implementing  a pitch  rate  control  system  for  the 
F-16  aircraft  clearly  indicate  that  digital  simulations  can 
provide  a very  cost  effective  approach  for  designing,  devel- 
oping, and  optimizing  advanced  aircraft  weapon  delivery  sys- 
tems. Used  in  conjunction  with  manned  simulation  efforts, 
the  digital  simulation  approach  can  provide  invaluable  in- 
formation. It  is  therefore  my  final  recommendation  that  both 
the  EASY  and  TAWDS  programming  techniques  be  considered  to 
aid  in  the  design  of  fighter  aircraft  handling  qualitites 
specifications  as  sot  forth  by  MIL-F-8785B. 
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Appendix  A 

EASY  Model  Program  Description 
and  Computer  Generated  Schematic 

Included  in  Appendix  A is  a description  of  the  EASY 
model  program  developed  to  simulate  the  F-16  control  system, 
pilot  mode],  and  aircraft  equations  of  motion. 

The  model  description  is  shown  on  page  95.  The  block 
diagram  location  of  each  standard  component  of  the  model 
description  is  specified  along  with  the  EASY  name  of  each 
component.  Also  shown  are  the  inputs  of  each  component  to 
specify  the  system  interconnections. 

The  computer  generated  block  diagram  is  shown  on  pages 
96  through  101.  The  block  diagram  verifies  the  component 
locations  and  interconnections  of  the  system  as  specified 
by  the  model  description  commands. 
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Appendix  B 

Baseline  Data  for  the  F-16  Aircraft  Simulation 

A portion  of  the  computer  output  of  Run  #43  of  the 
Griffin  program  is  shown  as  Appendix  B.  The  stability  axis 
derivatives  were  used  as  the  test  case  for  the  F-16  aircraft 
. simulation. 

Units  of  the  stability  derivatives  are  given  as  per 

degree.  However,  it  should  be  noted  that  this  applies  only 

. 

to  the  dynamic  derivatives  that  are  functions  of  primary  con- 
trol surface  deflections.  Static  derivatives  were  determined 
to  be  given  in  radian  measure. 

Roots  of  the  longitudinal  characteristic  equation  are 
shown  along  with  the  elevator  numerator  characteristics  that 
were  used  to  validate  the  dynamics  of  the  F-16  aircraft  model. 
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Appendix  C 

TAWDS  Programming  Details 

Included  as  Appendix  C is  the  Terminal  Aerial  Weapon 
Delivery  Simulation  program  data,  program  output,  and  addi- 
tional program  statements. 

The  extensive  data  list  required  for  the  program  opera- 
tion is  shown  on  pages  106  through  114.  Pages  106  through 
108  lists  the  input  data  required.  The  stability  derivative 
requirements  for  the  non-linear,  six-degree-of-f reedotn  simu- 
lation are  shown  on  pages  109  through  114. 

The  TAWDS  output  of  encounter  8 for  the  present  F-16 
flight  control  system  is  shown  on  pages  115  through  129. 
Pages  130  through  145  show  the  program  output  of  the  same 
encounter  with  the  pitch  rate  control  configuration. 

The  Appendix  is  concluded  with  a listing  of  the  pro- 
gram statements  required  to  adapt  the  generic  TAWDS  program 
for  an  F-16  aircraft  simulation,  pages  146  through  154.  The 
programmed  target  rates  are  shown  on  page  146  along  with  the 
target  maneuvers  of  page  147.  The  control  statements  added 
to  program  the  F-16  flight  control  system  begin  on  page  147. 
Included  are  the  logic  statements  for  the  gain  scheduled 
parameters  of  the  flight  control  system  that  are  functions 
of  static  and  dynamic  pressure. 
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no  rpt-MAT  ( I HO , TOM  ERROR TK.CNT  rx;F.iTnS  50CJ) 

• n*  | i i i 1 ■ | >j  1 . i ; 

• . ( ,i  ; OH 

•PH  m r,-inT . io  * 

RR‘-'l  02  ( 0)  = OHIO 
•OIL  FTC  PR  IN".  .'(• 

. ,?rrt.  1 . 1 10.1 ,5F8.r,F(>.0) 

• n^L  FT  : rtrifj-  .27; 

' 1 < 1 X,‘  H .0, 1 X,  Fa  .1 , 5F  8 .3,F8. 1, 7C8. 2,?rr,.2> 

• a i k'  r i i . 5ht> 

r>APftM(  6 0)  HHO 

• 0 ° 1 1.  T I 1.6*. 

m mt. - -ekeitl* 1000. 

•T)  n IL  T 1 1 . b5 

DTMT  -FRT  UL“  10P0  . 

•O'LFTr  T N l ) A r 1 .20 

.,  ( 5RFLTLC  , nilH0(  13)1 

•OFLt'L  TNPAC1.27 

.,  < F.PTRTLC  , niM3(  15)) 

• orL fte  tappet. sc 

. ai ,Tnco/-i.5,a.  o,  i.4, 2. a, 5, r.o  ,a.5,s.a,in.  o / 

• 0rL  STr  T ArCOrT  .S3 

. ,01 1'vnv-l  00.  .-80  . . -6C  . • -55  . , 0.  . 30  . , 55.  , 60  . ,80  . , 1 00.  / 

•Dim  TAKC.rT.n5 

. , v t*‘t  'i  /r’oip.0t0.o,o.c»o.o,n.o»o.0|0.o,o.o,o.o,o.o/ 

••'rLr:‘  TA^r.tT.n.; 

• n • L 1 F TCHMnR.<M  , 

. 5 5*; i 0 '3/3  , 

•OCLfTC  iCNH'-r.ni 

. 5^5 to  ’325, 

•oil  nr  iCNVF.nr 

. 55 5 P '372 » 

• 0 L E T r T C N f ! • 9 . 0 3 

. 55550 3372, 

•mnr  icnmcr.s5 

. Sm5.  C'3‘12, 

•OFETTC  TCNM!'P.)H 

. 555.01391, 

;priHr  iCNf'P.sa 

. 555-0  330.:, 

•OcLrK.:  TCNMrR.n’ 

. r>5  • 5 0 3391/ 

• n fit rr c icn*' r.icf. 

«0«, 8., IP. ,11. ,13., 13., 15. ,1"..  / 

• OFL  ( TF  ICNH'-i  . ( t,. 


f 


ffltl  PA 01  18  BIST  QUALITY  PRACTICABLE 
WfY  FURHISHKI)  IOW>C  . 


*0ELETF  ICNMPR.146 

ITMX(1,1>=E60  $ ITMX(1,2)=410  $ ITMX ( 1,4) =410  $ ITMX(1,F)=410 

•OFLFTF  I CNMr  R . 1 47 

ITMX ( 2* 1 > =37  0 $ ITMX(2,2>  =410  $ ITMX  (2,*>)  =310  $ IT  MX  ( 2 , 7 ) = 4 1 0 

•OFLFTE  TCNH<^.14M 

TTMX(3»1)=370  $ ITMX(3,2) =410  S ITMX ( 3,7) =610  S ITMX ( 3 , 8 ) = 41 0 

•OFLFTF  IClJMrR.14‘T 

ITMX(4,l)=o60  E ITMX(4,2) =410  $ I T MX  ( h , 4 ) * 6 1 3 S ITMX  (4 , 5)  =41  0 

•OCL^Tr  ICN^o.i.jj 

!TMX(?,l)=f70  ? ITMX(5,C)  =410  T ITMX  (5,4)  =P>10  ? ITMX  (5,5)  = 410 

* on.  r T r i 

lTM(f.l)  3;'<  ) t ITMX(6,2)=410  $ ITMX  (6,4)  = 310  $ I T MX  ( 6 , 5 ) = 41 0 

* DEL  E~F  I CN“r'R  • 1 S? 

ITMX(7,1)=380  I ITMX (7, 2) =4 10  $ I T MX ( 7 , 3 ) = 3 1 0 f ITMX ( 7 , 4 ) =4 1 0 

* DrL  FTr  ICNMpR.163 

TTMX ( 4 f i ) =o50  f JTMX(8,2)=410  S I T MX < 8 , ? > =610  $ ITMX ( 8 , 4 ) =41 0 

•OFLFTF  I C N m r R .171 
21  I T H ( I ) = 4 1 0 
•OTLETr  fUTSIl.2,365 

SUBROUTINE  AUTSIl 

COMMPN  /cruu/  n'l'U530),  OUM1  (270),  DUM?(300> 

COMMON  / 1 0 8 T A / 1 OUM  1 I "■  0 0 ) , lOJ-3(?C> 

COMMON  /ItlTVAR/  ’’AFAH  ( 46  0 ) , TIME,  IMTOCX 
COMMON  /CINT/  T.HMftX 
REAL  MfCH 

OI  MENTION  ALPHA  ( 4 ) , H <■»)  ,P(4)  , O ( 4 > » R ( 4 > ,1  NODR  (23) 

DIMENSION 

1 CpM7(D,  CCN7(?>,  DVNT(R),  COOQC’),  CCOQ(?>,  nVQO<8) 

?,Cp“0Ml  (3)  , CC  COML  ( 3 ) , PV/COML  (1?>,  COSF(2),  CCSC(2),  nvcrm 

3 , CP.RF  (?)  , CCRF (2 ) , OVRF (?) ,CnYRv2> , PCYR(2>,  0VYR<8) 

4, CONV (?)  ,CCNY  (?)  , PVNYtd),  DVER2(3),  )VFRS<3>,  0VPHX(8),  OVFDS(8) 
EOHIVALENCE 

A ( GKH72  ,PUM(10r>>,(  GKFR  ,0 JM ( 1 1 7) ) , < GKN71  ,0UM(U9>) 

R,(  GKNt  1 , OUM  (12  0 ) , ( GKMFTH  ,0114(125)), ( GKNL2  ,DUM(1?S>) 

C,(  GK6  ,r)UM(14?)  ),  ( -.KF  ,0  JM  ( 1 5 2) ) , ( GKNZ?  ,OUM(106)> 

FO"ivtLfNCE 

a (rpN7(«.)  ,nuM(io«)  >,  (cnN7<i)  ,o.f*(  it  d)  , (coonm  ,m)M(U3>) 

n,{croo<t>  ,ouM(iir) ),  (rnsEd)  , oj  > i * ? 1 1 > , tccse  ( n ,oum(i23)) 

C,  <nCO"L  (1)  ,OU«(  1 4 3 > > , (CCCOHL  ti) ,1J«( IB*) ) 


FOCI  VAl  F NC E 

( PR 

, Pl)M(  3)) 

EOMIVAI  ( MCE 

( MACH 

, OUM (038) ) 

EO'II  VAl  FNCF 

( GKQ 

, 0(11(191)) 

FONT  VAl FNCF 

( GGNF 

, OIJ1(42f)) 

EOUI VAl FNCF 

( GGNA 

, PI  H ( ♦ 2 6 ) ) 

EOUI VAl F MCE 

( OA 

, OUM 1 ( 94>) 

EOtllVA!  FNCF 

( oahi 

, OUM 1(112)) 

FQUIVAl FNCE 

( PALO 

, OUM 1(113)) 

FOUIVA1  FNCF 

( OSLO 

, OUM 1(114)) 

EOOIVAIF  nCF 

( OSLL 

, OUM  1(110) 

EOUI VALFNCE 

( PSHI 

, OUM 1(116)) 
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copy  fuktushisd  to  udc 


, 0^1(117)) 
DUM1(121)>  , (RUDLO 
, ou-uasg)) 

, 0UH1(1F.2)> 


0UM1 (1?2) ) 


, 0im(1561> 
PUM(6?C) ) 


EQUIVALENCE  ( OSLO 

EOUIVAl  F MCE  ( PUPHI 
cQUIVAiMlCF  ( PHIO 
EQUIVALENCE  ( OCOM 
*,  (P'.!OCrM,niMi  (165)  ) 

EOUIVAl FMCE  ( OOS 
EOUIVAl ENCE  ( FP1, 

EQUI V At  ENCE 

A C CBR-d)  ,QUM2  , < CCnr(l)  , D JM2  { i)i.  2 ) ) , ( COYR(l)  ,nUM2<CpG)) 

«),(  Cf.vc  ( 1)  ,OU"2(GrB))  • ( CHNY(l)  ,DUM2  (062) ) , ( CCNY(l)  , DUM2 ( 064)  ) 

EQUIVALENCE  ( IRUNNO,  IOUM3(32>) 

EQUIVALENCE  (IFRNT,I0UM1(132)) 

CQUI V AL  c NCE 

1 < ION7  , I NOOR ( 1))  , ( IOQQ  , I NOOR ( 2)),  ( IOSE  ,INQOR(  3)) 

2,  ( IPCPML,  I NOOR  ( P>>  , ( IOYR  , INnOR(12))t  < TONY  ,INQOR<13)> 

EQUIVALENCE 

A <nvER?(l> ,PARAH(l?g)>  , (OVERS (1) AR4 M ( 1 37 > ) , ( OV°H X ( 1 ) , P AR AM ( 145 ) > 
'f»,  (CVEPP  ( 1)  tPARAM  (163)  ) 


, QV  E°2 ( 5 ) ) 
, OV  EPS  ( 5 ) 1 
, OVPUX ( 5 ) ) 
, OVFPS (5 ) ) 


EQUIVALENCE  ( ER2X 
EQUIVALENCE  ( EkSX 
EQUIVALENCE  < Er MX 
EQUIVALENCE  ( FPSX 
EQUIVALENCE  (H(l) 

♦ ,(P(1) 

* , (R ( 1 ) 

DATA  nnR  / 57.2957795131/ 
DATA  ERRAUT/  8MAUTCI1  / 

RETURN 


( ER2I  , 0 VFR2 ( 1 ) ) 
( ERSI  , QVERS(l)  ) 
( PPHI  , QVPHX(l)  ) 
( PRSI  , 0 VE°S ( 1 ) ) 
,"ARAH(77)  ) , (ALPHA  ( 1) ,PAPAH( 17) ) 
,»ARA*f  (25)  ) , (Q(l)  , •>ARAH(  33)  ) 

, PARAM (41)  ) 


ENTRY  AUTS.I2 

CALL  I) UFO  (129, AOUMXX, IOUMXX) 
CALL  miPO  (133,  AOUMXX,  IQU^XY) 
CALL  INUPQ  (137, AOUMXX, IOUMXX) 
TALL  I'"PQ  ( 161  , AOUMXX,  IP'JMXX) 
CALL  . I»  IIPQ  (146,  AOUMXY,  nJMXy) 


CALL  I'UPO  (149, AOUMXX, IOUMXX) 
CALL  I>UPD  (153, AOUMXX , IOUMXX) 
CALL  I f UFO  (157, AOUMXX, IOUMXX) 
ISAVTR  = 1 
QO  ICO  I = 1,  23 
100  I NOOR ( I ) = -1 
NSN'IT  = 0 


RETURN 
ENTRY  fUTS 

IE  ( IS7VTR.NE.1  ) GO  TO  110 
IE  (GKNL2.LT.1.0 ) GO  TO  2 
1 1 - GKNl  2 


QEL-HMAX'ELOAT(I I> 
QELT  = 0 . 0 
GO  TO  * 
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f 


2 nELT=lC000C0. 

U rONTINt'F 

IF  ( I R'.'NNO  • GT  . 1 ) GO  TO  6 

C RECONVERT  LIMITS  FROM  RAOIANS  TO  DEGREES  (CONVERTED  TO  RAOIANS  IN 
C SincotiTIN*  OAT  A I N) 

OSHI r Of  H T*  OOR 
OSLOs  Of  L 0*  DDR 
OSlU^O'LU^nnR 

r,siL  = rr  i L*Uf)R  , 

DA“I-iU  HI* DDR 

OALO=0/‘LC*OOR 

PIJOHl  = 'UCMI*DOR 
RUOLO=tUnLO*DOR 
0STAR  = rns»nr)R 

•C INITIAL  CONDITIONS  • 

•f.  IF  ( OS” A ° . GT  . OSH  I ) OSTAP=DSHI 
IF  (OSTAt,.LT  • DSL 0)  0STA8=0SL0 
RSC1=ALDHA(1)»00R 
CS>C1=P‘  01*C9NZ(1> 

X5=CEG! 

R5C2=0<1 >*009 
C5C2=RrC2*CR00(l) 

xa=c5o." 

AZ3  = GG»'A*GK»172*Xi.*0KF 

C OS  IS  COMPACT  PRESSURE  PS  IS  STATIC  PRESSURE 
o«5  =PA  R A m <«♦  SO > * 1 7 1 F . 0 * ( F 1«  . 7- 0 . 0 0 3S 55 * H ( 1 > ) 

08- PS*  ( ( <G.2*MACH**2U.O)**7.5>-1.0> 

IE  (OP. GT. 3333.)  GKNZ1=.G53 
IF  (08. LE. 7300.)  GKNTi=-.0002*08*.5B3 

if  (OP.Lr.eoa.)  r.Kiri=-. 30089*00*1. 2; 

• IF  (08.LF.23j.)  GKNZ1=1.Q 
»X«.sGKf71*X4*  .7 
AU^=(-1S.0  AX4  ♦ XF)*G<F8 
IF  (AU7.LT.0.D)  AU^O.O 
TP  ( AU3 .GT  .999.0)  AU3  = 999.0 
U3  = A73»AI»3 
R5C8=tr 

C5fe=C8CONL(l)*R50fl 

FIX  TO  SFT  PITCH  INPUT  TO  ZERO  PY  RIASING  TRIM 
X 3 =CS  C 8. 

X1=X3 

C5T3=C? Cfl  , 

RF03=0.0 

IF  (Cr"r  (1) .N-.O .0)  P507-C50?/CPSE(l) 

UJ-SEC7 

OCOHO=(U1-.A4)/O.318+7.2F 
XlTPIMr-X3 
C503*0.0 
8503*0.0 


f* 

!! 
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ui=c.e 

OCOMO=C.O 

xi=o.o 

Axi=G.r 

AX5=  (X* +AX4-30.4  >*r-KN73 
IF  (AX'. LT. 0.0)  AXB=C.C 
IF  (Axr  .GT.919.D)  AXP=R99.G 
C C.KNL1  JS  FUNCTION  OF  DYNAMIC  PRFSSJRE,  0" 

C GAIN  F/CTOR  CF  3.G  IS  TNCLUDFO  IN  FUNCTION 

IF  (OB.GT.330J.>  GKNL1=  3 . * ( . f 8 3) 

IF  (CP. IF. 3000.)  GKNL1=3.* ( - . 0 003* OR ♦ . 68 3) 

IF  (Q9.LF.800.)  GKNL1=3.* (-.00089*33*1 .25) 

IF  (00. IE. 330.)  GKNL1=3.0 

AX6=(AX1»AX5)*GKNL1 
9ATI0=rB/PS 

IF  (RATIO. LE. .53)  TRFAC=.5 
IF  (PAYI0.GT..53)  TRFAC=-1.19*RATI0+1.13 

IF  ( FAT.  10.  GT  .1.79)  TRF  AC  = -1 • 0 
X6=X5»-RFAC 
FPFI rPf  T A 9 
X2=Er 5T-X6-AX6 
F»?I=X? 

IF  (GK* .NE.9.0)  ER2I =X2/SK6 
A X?=  G . C 

FPHI=0.0  ... 

FPPI=C.C 
ISAVTF  = 0 

IF  (IRI’NNO.GT.l)  GO  TO  110 
CALL  LT  N FS ( 1 4 ) 

WRITE  * I FRNT  » 30 ) 

30  FopmaT  ( 1H1 ,1X,26HAUGMFNTATICN  VARIABLE  DUMP,//) 

WRITE  (TP=NT,  31)  R5C1  ,C5 C 1 . R5C 3 , CD? 3 , =>5  C , CE  0 ? , ”5  3 8 , C5 OB 

31  ror  (AT  «lH0.r-X,GI  PS01,TX,fi’l  C01.3y,6m  RG0?,7X,f.,M  rso’.rx, 

A 6H  Ff03,?Xf£.H  CI>.i3,7X,6H  P50T,7X,OH  C5 0 H , IX , oE 1 3 .6 , /) 

WRITE (IopnT, 33)  AZ  3 , AU3 , AX  1 , AX5 , AX  3 , X 5 , OST  A R , X 3 
33  FOPMAT  (7X,6HAZ3  ,7X,6HAJ3  ,7X,6HAX1  ,7X,6HAX5  , 

A 7X.6HAX6  . ,7X,6HXf  ,7X,6HDSTAB  ,7X,6HX3  , / , IX , 6E 13 . 6 , / ) 

11C  CONTINUE 

AS  OF  TME  DATA  15  AUGUST,  1975  , THIS  PROGRAM  USES 

36F  LOCATIONS  OP  THE  PAF AM  ARRAY.  TMr  RPEAKOOWN  IS 
PILTI1  - 53  LOCATIONS 

AUTSI1  - 188  LOCATIONS 

AIPFT1  - 1 ?3  LOCATIONS' 

WITHOUT  ALTERING  THE  ° A RAM  ARRAY  AN0  CONSrOUENT  PROGRAM 

LOGIC i AN  ADDITIONAL  14  INTEGRATION  VAPIA9LES  CAN  BE  ADDED. 
TH*SE  VARIABLES  CAN  BE  INTEGRATED  BY  EXECUTING  THF  CALL  TO 
EITHER  INTEG  PR  TRANFR.  IF  ONE  CHOOScS  TO  CALL  INTEG,  THEN 

THE  PAPAMETER  INTDEX  MUST  BE  INCREASED  9Y  8 POR  EACH 

ADDITIONAL  CALL  TO  INTEG. THAT  IS  INCORPORATED  IN  THE  PROGRAM. 
T HT  PApAMCTrR  INTIi-X  IS  INITIALISED  IN  THE  AIPTI1  SUD.RO’ITtUF. 
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C 

C« 


c 


L0*'GITU9INAL  CHANNEL  AUGMENTATION 


U5=ALPHA  (1)*0DP 
IF  (UC.GT«30.0)  U5=3C.0 
IF  (U5.LT.-5.0)  U5=-5. 0 


call 

TF  AHFR ( 2, 

COf'Z,  2, 

CCN7 , 

8,  9VNZ , 

ERRAUT , 

501,' 

U5» 

A 

T9U7, 

Y5,  r5  C t « 

0,0, 

ssoi. 

0. 

0» 

UUrrj  ( j > » ;(j' 

call 

tfahfm  ?, 

C!1PO,  2, 

CCQO, 

a,  3V0Q, 

ERRAUT, 

902, 

U4 , 

A 

1099 , 

X4 » 050  2, 

0.0, 

C502 , 

0, 

0) 

A??sGGf'A»GKN72  + Xi4*CKF 

09  IS  COMPfCf  PRESSURE  PS  IS  STATIC  PRESSURE 
PS=PARAM  (480)*1715.0M518  .7-0 . 0 0 3555*  H ( 1 ) ) 
09=PS*  ( ( < C.2*MACH**2*1.  0 ) ** 3 .5 ) - 1. 0 ) 

IF  (On.GT  .3000.)  GFfi?i=.C83 

IF  (OB.LE.30C0.)  GKNZ1=-. 0002*09*. 683 

IF  (QO.LE.80 0.)  GKN71=-. 00089*09*1.25 

IF  (OP.LE.29u.)  GFN71=1.0 

ax*,  = GKf'?  1*X4*  .7 

AU3=(-15.0  ♦ AXu  * X1')  *G<P9 


IF  (AU’.LT.O.u)  All  3=  0 . 0 
IF  (AU’.GT  .999.0)  AU3=9°9.G 
»3=AZ3*AU3 

CALL  TP AfIFR  ( 3,  C3COML,  3,  CCCOML,  12,  0VC0ML, 
A inCOML  , X3,  R5 P8,  0.0,  C50S,  0.  0) 

FPlsOCCN+OCOMO 

IF  (FOl.C-T.  7.25)  Ul=.318*(FPl-7.25)  ♦ • 44 

REMOVE  OE A 99 AND  9ETWEFN  -1.75  AND  1.75 
IF  (FP1.LE.  7.25)  II 1= . 06  0 69*  FP1 

IJ1= . 08*  (FOl-1.75) 

U1  = C.  ’ , 

Ul=.  06*  (PP1*1 .75) 

Ut-.’42*  (PP1*7.?3)  -.44 


IF  ( FP1  • LE  7 .25  ) 
IF  (FPl.LE,  1.75) 
IF  (FPi,tE.-1.75) 
IF  (FPi.Lc.-7 ,2") 


IF  (Ul.LT.— «. 0)  UI--4.0 


ERRAUT,  508,  U3, 


IF  (Ul.GT.8.0)  Ul=8.0 

CALL  TFANFRl  2,C9SE,  2,  CCSE,  8,  9 VSE,  ERRAUT,  503,  Ul, 
A IPSE  , XI,  P50  3,  0.0  , C50  3 , 0 .0) 

AX1  = X3-X1*X1To.IM 

AX5= (Xf *AXU-20.4)*C<N73 

IF  (AXr .LT.3 .0)  A X5  = C .0 

TF  (AX*  .GT. 199.0 ) Axr=9°9.0 

GKNL1  ?S  FUNCTION  OF  DYNAMIC  PRESSURE,  09 

GAIN  FfCTOR  OF  3.0  IS  INCLUDED  IN  FUNCTION 

IF  (OR. GT. 3000  .)  GKNL1= 3 . * ( . C 8 3) 

IF  (09. Lc. 3230.)  GKNLlrJ.* (-.0002*93*. 683) 

IE  (r  9.  LF . " i ; C.KNLls3.*(-. 00081*09*1. 25)  • 

IF  (OB.LE.78J.)  GKNL1=  3.0 

AX6*(AX1*AX5)*GKNL1 

RATIOa*'9/PS 

TF  (RATIO. LF. .53)  TRFAC«.«5 
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TP  (RATIO. GT, .53)  TPFAC=-1 . 19*RATIO*l . 13 

IP  (RA1IO.GT .1.79)  TRFAC=-1.0 

X6=X5*TRFAC 

X7=AX6»X2 

X8=X6»X7 

IP  (XH.GT.nSLU)  GO  TO  40  ’ 

IP  (X8.LT.PSLL)  GO  TO  42 
X9*3  • 0 

GO  TO  '.4  , ' 

40  X9rG‘'6MXB*0SLL) 

GO  TO  w4 

42  xn-GK6*  ixa»nsLU) 

44  U2  = AXf,  - X9-AX2 
F.R2X  = U2 

CALL  IfTEG<E02X,ER2I) 

X2=GK6*ER2I 

IP  (X2.GT.0SUD  GO  TO  46 
IF  (X 2 ,LT . nSLL)  GO  TO  48 
A X2- 0 . ( 

GO  TO  40 

46  1X2  = 2010  .* (X2  + PSLL) 

GO  TO  < 9 

48  AX?=20C 0 .* (X2+0SLU) 

49  rcnii1'? 

STAC0MrX6+X7 

RST9=2C  .0*  (STACOM-FRSI) 

IF  (FST9.GT.  60.0)  RSTn=  60.0 
IP  (RS*9.LT.-60.0)  RSTB=-60 • 0 
EPSX=RrTR 

CALL  I(!TEG(ERSX,ERSI) 

OST=ER?I  , 

IP  (DST.GT.OSHI)  D?T=OSHI 
IF  (OST .LT.OSLO)  DST  =DSLO 
OOR=OST/noR 

6 LATERAL  - PIRfCTIOHAL  A'lTOEHT ATION 

C- 

C Alirpnn  FHANtlEL  COMMANOS 
ca  = PHI<~  * 

IP  (FA. LT. -11.0)  PHICl=?^.0*FA+a38.0 
IF  (FA. GE. -11.0)  PHIC1=12. Q*cA  + 52.0 
IF  (FA.GF.  — 6.0)  PHIC1*  3. 33?333333*FA 
IF  ( F A ,GC.  6.0)  °HIC1=12.0*FA-?2.0 
IF  (FA.GT.  11.0)  PHICl=33.G*FA-?38.0 

CALL  TF  ANFF ( 2,  CORE,  2,  CCRE,  8,  DtfRE,  ERRAUT , 511,  PHIC1, 
• IORE  , PHTC2,  0.,  0.,  0.,  0.) 

PHIC3*C.12*<P<1>  *OPP-PHIC2) 

FPHA«2C . C* (PHIC3-FFHI) 

IF  (FPHA.LT, -80.0)  FPHAs-BO.O 
IF  (FPHA.GT.  80.0)  FPHA-  80.0 
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UOrl  ro  2D£  ^ 

FPHX=ErHA 

CALL  IfTEG(FPHX,FPHI) 
pn»=FPuI 

IF  (0(V  iGT.OAHI)  DPA=OAHI 
Tf  ir'.n.^v.ni  r>r  r.  -tplp 
f>6  : r-C'rt  / 

2 <unr>tr'  channel  ( onmaniis 
rupci=‘  upcom 

IP  (RUrCi.GE.  15.0  RJPC?='0 . ?16*RLn:i-64.75 
IF  (RUrCl.LT.  15.0)  R'J"C2  = 0.0 
IF  (RUPC1.LE.-15.0)  R'JPC2=0.316*RU3:i*64.76 
YAWERsf (1) *DPR-X5*P(1) 

Rl=R(l) 

01=0(1) 


CALL 

T f ANFP.  ( 2 

,CBNY, 

2,  CCNY, 

8 , 3VNY  , 

ERRAUT  , 513  , YAWFB , 

• 

IPNY 

, RFP1 

i C • i 0 . f 

0.,  0.  ) 

• CALL 

TFANFF  ( 

?i C"vp 

, 2 , CCYP, 

8 , OYYR  , 

FRRAUT,  512,  RFB1 , 

* 

IPYR 

, RFB2 

, o . , c . , 

0 • f 0 • ) 

F6  = PF“?*C.f)*Gr.HS 


IF  (RA’IO.GT.3.3)  AKE9=F8 

IF  (PATIO. LF. 3 .3)  AKF5=-F8*  ( 0 . 435*RA T 10- 0 . 30 5) 
IF  (RA*IO.LE.2.0)  AKF3=0.5*F6 
RinOC  = f KF9-^onC2 
FPSA=  2„ . C* (RtnnC-FPSI) 

IF  (FP^A.GT.  120.0)  FPFA=  120.0 
IF  (FPFA.LT. -120 .0)  FPPA=-120.0 
FPFX=FFSA 


CALL  r»'TEG<Ff'FX,  FPFI) 

.luf-rr  I 

IF  ( OF « ,i.T  .R:n‘<l  I CPRt  Pl'PKI 
IF  (0RC  . L I . RULILO  ) DP.R=  RUOLO 
nREDPP/PPR 


IF  (IRt'MNO.GT  .1)  GO  TO  80 
IF  (PFLT.GT.T)  GO  to  TO 
60  IF  (PEI T.GT.C.O)  GC  TO  65 
WRITE (JF=NT, 51) 

51  FORMAT ( 1*1 ,?Sm1AUGMENtATI0N  TIMr  HI S T OR V , // , 5X , 6WT I ME 


A t7X,£-PGGWA  , 7 X , 6 HU4  ,7X,6WU5 

R 7X,fcHSTAP0M,  7X  ,C.M£CSX  ,7X,6HER5I 
WRITF ( ]PPNT,53) 

53  FORMAT  (19X,6H°HIC1  ,7X,6H=>1 

A , 7 X. f ►•fphi  ,7X,oMR1  ,7X,6HYAWF9 
WRITE ( tdcNt,57) 


,7X,6HU1 


i 7X, 5HAX6  ,7X,6HX6  , 

> 

» 7X , 6HPHIC2  ,7X,6HPHIC3  , 
.rx.SHRFOl  , 7X,  4HRFP2  ) 


57  fqpmaT  (18X,6HF8  ,7X,6HA<fj  , 7X , 6MRUOCOM , 7X , 6HRU0C 2 , 

A 7 X, 6UX5  ,7X,6HRU00C  ,7X,6HF3SA  ,7X,5HPHSI  ,7X,6HGGNS  ) 
WRITE (TPPNT, 59) 

59  FORMAT  (18X,6HOS  ,7X,6H0:  ,//> 

65  OELT=OrLT ♦OFl 

70  WOTTf  (tppnT,71)  T,  111 , GG NA , 04 , U5» AX6, X6 , STACOM, RSTB , OST 
M FORMAT (IX, 10 E 1 3 . 6) 
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WRITE  ( JPPNT  , M)  PHICl,Pl,«>Hir2,PHi:3,DDAtPl,  YAWFB,PFB1,RFR2 
73  FOPMAT(H.X,9E13.6) 

WRITE (JPRNT, 73)  F8 » A Kc8 , RUOCOM ,RUD32  ,X5»RU0DC» FPSA , DRR, GGNS 
WRITF (IPRNT,77)  PS, OR 
77  POFMAT U-.X,2£13.6) 

SO  CONTINUE  ' 

RETURN 

END 
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Appendix  D 

EASY  Analysis  Program  Data 

The  program  commands  of  the  EASY  Analysis  program  are 
included  in  Appendix  D. 

Table  parameters  specify  the  independent  and  dependent 
variables  for  the  table  look-up  gain  parameters  of  the  system 
model.  The  parameter  values  which  satisfy  the  input  re- 
quirements of  the  standard  components  are  listed  following 
the  tabular  entries. 

The  longitudinal  axis  stability  derivatives  for  the  F-16 
aircraft  at  the  selected  flight  condition  of  .8  Mach  and 
20,000  feet  are  shown  on  page  157.  In  addition,  page  157 
shows  the  commands  necessary  for  generating  a steady  state 
system  solution  and  establishing  the  aircraft  trim  condition. 

The  frequency  domain  analysis  is  completed  with  program 
commands  to  establish  a pseudo  tracking  task.  Program  com- 
mands for  both  closed  loop  and  open  loop  analysis  are  com- 
pleted on  pages  157  and  158. 

The  EASY  Analysis  program  data  list  is  concluded  with 
program  commands  to  generate  a closed  loop  system  time 
response  to  a reference  step  input. 


THIS 

FROM 


page  is  best  quality  practicable 

COPY  FURNISHED  TO  


T A TL  r»  FT."  "UEi  , - 
-1~  .'.v,  -"  .25 ,7  ,2S*-3  . 

- <i . , - . 4 4 , • *s  4 , 1 0 • 3 •') 

TA3Lr,  FT  AFUFE,  w 

“1  • , “i  i , — 4 . , —4  • 

T/''>Lr»  "TAFIIE?,^ 

0 i » ? ^ J i , 5 j w .,7'«C««,l>CG&* 

• 7 j » i , •37**l,.u541,.u531 
TA^L  ",  5TA  C"E‘  , r- 

d . , 2 o :• . , ° o : . , 3 j c j . , ft  c o c . 

3.  ,1 . , 1.  *"->“,  .249, . G^dl 
T AIL  r , FT  A “UF 1 ,4 
C . , . c7, 1 .79  . C 

• ,™2*G 
TAX^GTAFUFF,- 
•50. , -25 . ,25  . , 50 . 

-25.  ,0.  ,C. 

DA->V'?TEC  VALUES 


AN 

cUFi=-l, 

Af 

C|JF2  = -1 , AN 

FNE3=-1,4N  nJE+=-l 

AN 

~UF5  = -1, 

5'JE6=1 

Cl 

•i  a =•  i = i . 

Cl 

54-2=. 5, 

C? 

SAF2  = S . , C4 

S AE2  = . 5 

C3 

,56 

SAE3=-3  . 

▼0 

LG" 1=1 0 , 

,rG 

L 5i  r.  1 = 1 u • 

G A IL  FE2=  3 . , 

7r 

LEE2=Q.,C'0 

L EE?  = 1 . 

C2 

NCrl  = J . , 

C7 

UCE1=1 . , C“ 

MCE 1=- . , X3  MCEl=-?G  . *♦ 

C 2 

MC“2= 1 . , 

c: 

MC  E 2 = 1 . , C4 

MCE 2 = 0 • , X 3 ‘*1rE2  = - 15  . 

CF 

3 4 E *♦  = C . 

Cl 

■1Cr3=  , 1 6 1 , C 

7 yCE7=.16~ 

, C3  MCE  3 = 0 . 5 , C-+  MCE?  = 

Cl 

3 A " = .3, 

pr 

5,".E 3=3  . 

GATLEF"=3.C 

,7c 

LEE3=‘'«,pC 

L EF3=1 2 . 

Cl 

cr 

MCE4=1. ,03 

M5.i->  = -l . , C4  MCZ,=  0. 

c? 

ici:=-3. 

*53 

MCZ3=u.  ,Ci- 

MCC5=0 • , Y3  ^CE5  = 0 . 

GKHT-l  = r.  , 

GXIT- 1=200  0 . , 

AMAI7E1=2E. , AMIITEi=- 

n o 

U w. 

S1  4 E 2 = 0 . 

C 3 

‘1C“c=0  . , 

C- 

HCE3=0  . , X3 

MCEG  =C  . 

Cl 

M A F 3 = 1 • 

Cl 

MA^4=1  . 

Z1 

Tf-l=C., 

7 £ 

TFEl=27w-,«  , 

pi  TFri=-,2.a,°]  Trri  = 

Cl 

MAP3=-1. 

Cl 

5AF1=?3 . 

,57 

SAr 1=3 . , C4 

FAFl=20.,C6  5\rl=-3. 

7 C 

LGr 1= 1 . , 

F C 

LGF 1=0. 

C 3 

5Ar2  = ?‘>. 

» 5 6 

5?F2=-25. 

Cl 

MAE2=.0c 

U 2 

32 

156 


THIS  PAGO''  IS  PEST  QV Ml  TV  mCTICA^I 
COPY  tUvNlSH^J  TV'  UDC 


GAiLEnL=?4p .t u:  ?o  lfpl=o.c  pc  lfdl=20.c 

Z* ZL APL  = ?. 1M 5 7C  l A pL= . C 

ci  °f pl  = i . or -:•=>, c?  sapl=i.3,C'  s'\ri=o.:,:‘*  sapi  =1 . cE-ao 
C5  = cp  oapl=:.c 

71  TFPL  = C.C  7:  T FPL  = 1 . C PI  TppL=1.2  PS  TFPL  = t.C 

Cl  N0PL=1.C  C2  MCPL  =1.0  n2  MCPL=1.:  C'*  MCpL=0.0 


7 G LGpl=.?r 

PO  L GPL =0.0 

Cl  1 ApL  = -l , t 

O 2 MApL  = r:.5 

TX  S 0=0,, VO  r 

D=0.,T7  S^=C. 

TXX‘in=  1G  C7  . F 

rYYE0=4P9G-'>. 

17790=5^ 

7-0  . 

ix7C:o=  ipp  , 

I Dl  l V=3  . 

VS  pq.L 

1LSAV=2. 1 07° 

S AV  = ? 

3 3 • 

XC  LO=-.J2r;C 

XA  LO=-.‘“2-j1 

XU  LO=-.0 

7-+& 

Xr>ELP=.Cr2'> 

70  L 0= - . 1 3 

7 A LO=--,.»l'>p 

7ADL0=.  i-i 

] o 

70  L0=-2.5Hf>5 

711  LO=-.l?lc 

7.0910=  -.-.9  si 

MG  L0=-.I 

132 

MA  L L 0=  . J ° 

M A 0L 0=- .p  3D C 

MO  L0=-2.31i7 

MU  L0=- . 0 

145 

MOELO=- . G65P 

MA  1 L 0=  ••  P r ,5 

C 10=11.32 

XP1L0=3. C,F?1 

LC=J,TY1L0=C 

Cl  1C7 1 = . 0 3 1 1 

,02  MC7l=-.G311,C? 

MC71=G. ,C4 

MC  7 1 = 

- 1 * X 3 MC  71  = J . 

Cl  N A r T = - 1 t >r 

2 maet=?o jo r . 

7 C LPET=.1,F. 

LGET  = L • 

Cl  M A E N = - 1 . , ? M A £ M = 3 2 p t r » 7 C LG 

EN=10.,P3  L 

G E N = 0 . 

7C  lG-?=6.  ?,F 

0 LOE 2=3.3 

TNT  CONTROL? 

V =59=3.,°  f 

D = 0.,R  S0=C.  »7OLP0  = O . , YAWS'*  = 

3 . 

F p 7 0 R COMPfltS:  H S 3=  . 6 » W 0O=.f6,Q  $i  = . IE-03 

X2  LGEN=4.,\?  i.GFT=  .2F-C5,  X2  LG?2s.P"-;i»X2  l 3r  1 = .4r-C  2 
X?  TFP1=.'01,X2  LGF1  = .CCl,PITGC)=.4E-j?,UTSP  = 20. 

XT  LEE3=  .?F-it,X2  ITEl = . IT  4, X I TFFlr.i 
XI  LErL=.;,X2  L A PL  = . o C r- , X I ’ FPL  = . C Cr-  , *2  TFPL=.00r> 

X7  L G r L = • l 0 2 » X I ICE2=.4F-C5 
INITIAL  CONDITIONS 
ft  L T-^ °=  ^ a C- 3 f , » SP=3  29.5 

PRTNT  C ONT P Cl  = 3 
PLTr  ON 
PRINTrR  FLOTf- 
IN"  CONTROL = y?  LGpL  =0 
5 T E A OY  STATF 
XIO-X 

INF  CO^'T^OL  = X?  LGr N=0  > 72  L G F T = 0 
INT  COHTPOL=  X 2 LGPL=1 
LlNEAp  ANAlYcIc 

TITLF=thE7A/'MFTA-?EP  CLOc£P  LOOP  NORMAL  ACCEL 
TF  I M°UT  = C 2 MApL 
TF  OUT  PUI  = PTTSC 


157 


— 9 

sm  mt 

flflpy  |\^USUEP  15  i=M 


TF  '1AM UAL  SCALES 
FREO  MTN=.CG! 

FPEO  M \ X = 1 v.  C . 

300E  , T ’AVSFEF  FUNCTION 

TITLF='KFTA/*HETA=>FF  OPEN  LOOP  MORTAL  ACCEL 

DARAMrTEc  y/fiiUES  = Cl  MAPL=C.O 

TF  INPUT-  C?  ^ A PL 

TF  OU~PU*=  PJTfD 

TF  MAMMAL  SCALES 

FRCO  MTN=. C 1 

frfo  max=i:o. 

ROOF , T p A N S P F c FtJNCTIOM 
PARA  METE®  VALUE3=  Cl  MAPL=-1.0 
TTMC=.  C-5 
tm  \X  = R ‘ . *. 

0UTRAT“  = 1 . i 
PPATF=2i 
I MT  MOOE=E 

TTTLF  = CL03ED  LOOP  STEP  RESPONSE  • 
niSPLAYl 


oi-sn 

VS 

TIME 

o >n 

VS 

TIME 

A LT  SO 

VS 

TIME 

X 2 L OF i 

VS 

TIME 

SIMULATE 

! 


r 


Vita 

Michael  Marchand  was  born  in  Gonzales,  Louisiana  on 
April  6,  1948.  He  attended  Gonzales  High  School  there  and 
graduated  as  v?l edictorian  of  his  class  in  1966.  He  entered 
undergraduate  studies  at  Louisiana  State  University  and 
received  his  B.S.  degree  in  Electrical  Engineering  and  his 
ROTC  Air  Force  commission  in  1971.  Later  that  year,  he 
began  active  duty  in  the  Air  Force  as  an  Undergraduate  - 
Pilot  Training  student  at  Laughlin  AFB,  Texas.  After  re- 
ceiving his  wings,  he  attended  Pilot  Instructor  Training 
at  Randolph  AFB  and  then  spent  the  next  two  years  at  Laugh- 
lin as  an  instructor  in  the  T-37  aircraft.  From  1974-76,  he 
enjoyed  a tour  at  Mather  AFB  as  an  instructor  pilot  in 
Undergraduate  Navigator  Training,  being  a member  of  the  first 
T-37  squadron.  While  at  Mather,  he  accepted  the  additional 
duty  of  Functional  Check  Flight  pilot  for  the  T-37  mainten- 
ance squadron.  In  1976,  he  entered  the  Air  Force  Institute 
of  Technology  at  Wright-Patterson  AFB  to  attain  a Masters 
degree  in  Electrical  Engineering,  specializing  in  aircraft 
guidance  and  control.  He  is  married  and  has  two  children. 

Permanent  address:  Route  4 Box  245-M 

Gonzales,  Louisiana  70737 


159 


_ UNCLASS 1 FT ED 

5f  fUH|1  > Cl  A-  m»  l>  ■'  ' i > ' N «'•  i mis  I'AGE  /’*f#  / 


REPORT  DOCUMENTATION  PAGE 

NI'AD  INSTKUC  HONS 

Ml  I OKK  COMPI.I  riNf',  I OKM 

\ Ml  POM  1 NNMIll  M . i.o V 1 Av  v.  t bMON  NO 

V 

AF1T/CGO/FK/77-7 

» RF.CII'ILN  1 S f*  AT  AliH.  NUMtll  M 

1.  1 1 T LI  i «n.»  NwMlfl#) 

PITCH  KATE  FLIGirr  CONTROL 

FOR  TUF.  F-le  AIRCRAFT  TO 

S TYI'I  O f Ml  PORI  A PI  HIOO  ‘ ON.  1 Ml 

MS  Thesis 

IMPROVE  AIR-TO-AIR  COMBAT 

6 PCMI  OMMINl.  OMG.  Ml  I’OMI  NUM11LK 

7.  All  1 MOM,  *> 

Miehael  A.  Mai  chand 

Capt  USAF 

B CON  1 M ACT  OR  GRANT  NUMUI  M(*J 

>»  n HI  v'RMINk  ORGANIZATION  name  and  aoumess 

Ail  For eo  Institute  ol  Technology ( AF1T-EN ) 
Wr  .ioh  t -Pa  t t oi  son  AFB,  Ohio  -U'T3o 

10.  PMOGR  AM  ritMENT.Plh'JFCT.  T A‘>K 
AREA  A «HMK  UNIT  NUMMl.RS 

Project  7071 -00-12 

II.  a'MTROLUNw  Oil  ul  N AMI  AND  IL'I'MESS 

Air  f«m  vi  j.iia  Pyn.imics  r.ab"i  atory  vou; 

1.’  REPOR1  DAI  1 

December,  1977 

Mriahi  v.ii  t >.■'  avi>  i.'ii  4!'-l  ! ' 

11.  NIIMIU  M Of'  PAGES 

173 

14.  MONIlCKINu  AGENCN  N AM l A ADDMt.S''|d  dlllerent  from  Controliu i«;  iUlice) 

iv  SECURITY  CLASS,  (of  fbln  re, v?rP 

Unci  assi.f  i cd 

16*  Ol  a ASMnCAIION  I'OWNoU  APINv 
SC  HI  IHH  1 

16  CM  ST  RIHIl  1 ION  SI  Ml  Ml  NT  (of  thin  KV|».*rf ) 

Approved  lot  public  release;  distribution  i 

uni  i tn i t od 

1 7 . DIST HHUITION  STAU.Ml.Nl  (of  fh#  •/■»fmt*f  #nf#r»«f  in  Utoch  JO,  II  ditlerml  from  Neporf) 


16.  SU»’»’L«  MI.N  T ARY  NOT  LS 


Approved  vf  oi\oubl\Lc  release;  1AW  AFR  190-17 

^ V r,  'V  \ 

Jerral  F.  Guess';  Ca^fai  n,  USAF 


Director  of  Tn foi motion 

19  KLY  WOMOS  (Certf /mi#  on  r#v#m#  aid*  It  n#o#*.«*rv  end  Identify  by  block  number) 


i o,  Afft TRACT  (tontlnne  «vi  r#v#ra#  #M#  naceaaary  **1*1  Identity  by  block  number) 


■ ^ r — 

Digital  simulations  wore  developed  to  implement  a pitch  rate 
control  system  for  the  F-lb  aircraft  engaged  in  aerial  gunnery. 

I . r t,  t ho  EASY  Modelling  and  Analysis  Program  by  Boeing  Computer 
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performance.  The  Terminal  Aerial  Weapon  Delivery  Simulation 
(TAWDS)  program  by  McDonnell  Douglas  Corporation  was  adapted  for 
the  F-16  aircraft.  A non-linear,  s ix-degree-of-f reedom  aircraft 
model,  multi-axis  flight  control  system,  and  multi-axis  pilot 
model  were  developed  to  demonstrate  target  tracking  capabilities. 
Sight  different  air-to-air  scenarios  were  developed  to  simulate 
evasive  encounters  with  an  F-4  target  aircraft.  Time  history 
target  tracking  errors  indicated  the  improved  tracking  performance 
of  the  proposed  pitch  rate  flight  control  configuration  over  the 
present  normal  acceleration  configuration  of  the  F-16  aircraft. 


